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Abstract. We begin a study of a multi-parameter family of Cauchy initial- value problems for the modified 
nonlinear Schrodinger equation, analyzing the solution in the semiclassical limit. We use the inverse scat- 
tering transform for this equation, along with the steepest descent method of Deift and Zhou. The initial 
conditions are selected both to allow all relevant scattering data to be calculated without approximation 
y—( and also to place the governing equation in a transsonic state in which the quantum fluid dynamical system 

formally approximating it is of hyperbolic type for some x and of elliptic type for other x. Our main result 
is a global approximation theorem valid in a maximal space-time region connected to the elliptic part of the 
initial data. 

< 

^ 1. Introduction 

i— | Let e and a be positive parameters. The modified nonlinear Schrodinger (MNLS) equation 

^ (i) ief t +^0 + l^ + ^(l^)=o. *em, t>o 

is a completely integrable generalization of the focusing nonlinear Schrodinger equation, to which the MNLS 
equation reduces upon setting a ~ 0. The perturbation term proportional to a models the effect of nonlinear 
dispersion and is one of a suite of three terms arising in small-amplitude perturbation theory one order beyond 
I the focusing nonlinear Schrodinger equation in the analysis of short pulse propagation in weakly nonlinear 

J> optical fibers (TTJ (T2] [13]. Although the MNLS equation appears to be a perturbation of the focusing 

t-H nonlinear Schrodinger equation, there is a sense in which it can also be considered to be a perturbation 

of the defocusing nonlinear Schrodinger equation. This is related to broken Galilean symmetry introduced 
for a^O, which also shows that the MNLS equation can be reduced to the so-called derivative nonlinear 
Schrodinger equation of Kaup and Newell [15] . The complete integrability of the MNLS equation is connected 
with its representation as the compatibility condition of a Lax pair based on the WKI spectral problem of 
Wadati, Konno, and fchikawa |22j . The existence of a Lax pair representation makes it possible to analyze 
quite general solutions of the MNLS equation with remarkable accuracy; as an example, see the long-time 
asymptotic analysis of Kitaev and Vartanian [TO] [T7] . 
. £^ The Cauchy initial- value problem for ([!]) is to find a solution 4> subject to the initial condition: 

(2) (t){x,0) = A (x)e lSo ^/ e 1 xeR. 

Here, A (-) and S {-) are real-valued amplitude and phase functions respectively. Under suitable conditions 
on these two functions, the Cauchy initial-value problem can be studied with the help of an inverse-scattering 
transform derived from Lax pair representation of the MNLS equation ([!]). 

Our main interest is in the asymptotic behavior of 4> = 4> e ( x j t) in the semiclassical limit where e \. with 
a > and the functions A (-) and Sq(-) held fixed. This limit is obviously very singular, but the problem 
can be recast in a form that appears somewhat more tractable by introducing Madelung's fluid dynamical 
variables [15] : 

p e (x,t) :— \4> e (x, t)\ 2 , (a quantum fluid density) 

(3) f d I 

uAx,t) := < — \og(<fiJx, t)) > , (a quantum fluid velocity). 
I ox J 
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Here of course \og((j) e (x,t)) is denned (assuming p e (x,t) is non-vanishing) by continuation to be a diffcrcn- 
tiable function of x. The initial values of these fields are independent of e, since according to |2]), 

(4) Pe(x,0) = p (x) := A {x) 2 and u e (x, 0) = u (x) := S' (x). 

Most of what we will discuss in this paper has to do with quite special choices for the functions Aq(-) and 



So(-). However, to ensure the validity of certain general asymptotic formulae for scattering data (see (46) 



and (61)), we assume that the functions po(-) and u' (-) are Schwartz-class real-analytic functions of x £ 
and that po(-) is nowhere vanishing. In particular, these assumptions imply that Uq(x) — >■ u± as x — ¥ ±oo. 
Then 

(5) Sq(x) — S Q (Q) + / u (y)dy S Q (x) = u±x + S± + o(l), x ->• ±oo, 

Jo 

where 

p+oo />0 

(6) S+ := 5 (0) + / [u (y) - u+] and S_ := S„(0) - / [u Q (y) - «_] dy. 

In terms of the Madelung fields defined by ([3| , the MNLS equation ([!]) can be rewritten without approx- 
imation in the form of a coupled system: 

dp e i d ( _ , 3 2 \ du e d fl 2 \ 1 2 <9F[p e ] 



where -F[p e ] is given by the differential rational expression 

1 9 2 p, / 1 ^Pe 



(8) F[Pe] -2 Pe ^ \2p e dx. 

The only explicit dependence on e in either the equations of motion ([7| or the initial conditions Q lies 
in the dispersive correction term on the right-hand side of the equation governing u t in ([7]), and it seems 
reasonable to try to neglect that term and solve the corresponding quasilinear system of local conservation 
laws, the dispersionless MNLS system 

(9) | + £(^+^ 2 ) =0 and d it + l{\ u2 - p+apu )=^ 



with the e- independent initial data p(x, 0) = poix) and u{x, 0) = uo(x) given byQ. Although this seems 
attractive, proving that such a procedure yields an accurate approximation when e <C 1 is not at all straight- 
forward for several reasons. Perhaps the greatest obstruction lies in the fact that the dispersionless MNLS 
system (|9| is not strictly hyperbolic, leading to a certain ill-posedness of the e-independent Cauchy initial 
value problem purported to approximate the true dynamics. This means that without strong assumptions on 
the initial data, there may exist no corresponding solution of ^ at all. Even if there is a solution, lineariza- 
tion about a constant state reveals unbounded exponential growth rates corresponding to modes that are 
likely to be seeded by the dispersive correction term Q ; it is not obvious at all in such a situation whether 
the solution for small nonzero e should resemble for t strictly positive that obtained by simply setting e = 0. 

The characteristic velocities of the dispersionless MNLS system ^ satisfy a quadratic equation with 
real coefficients depending on the local values of the Madelung variables p and u. The discriminant of this 
quadratic is proportional via a positive factor to 

(10) Q := a 2 p + au- 1. 

At any given t > the fields p(-. t) and u(-,t) can cause the system to be in any one of three different states 
(with terminology borrowed from the language of stationary flows in gas dynamics): 

• Globally supersonic. In this case the dispersionless MNLS system is strictly hyperbolic, corresponding 
to Q > and hence real distinct characteristic velocities, for all ieK. 

• Globally subsonic. In this case, the dispersionless MNLS system is elliptic, corresponding to Q < 
and hence distinct complex-conjugate characteristic velocities, for all 

• Transsonic. In this case, the dispersionless MNLS system is hyperbolic for some values of x and 
elliptic for others, and Q changes sign as a function of x. 
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While the correspondence is not the obvious one, it is shown in [TO] that for initial data Q that makes 
the dispersionless MNLS system ([9| globally supersonic (respectively subsonic) at t = 0, the semiclassical 
analysis of the MNLS Cauchy problem can essentially be reduced to that of the defocusing (respectively 
focusing) cubic nonlinear Schrodinger equation. This result generalizes to arbitrary genus a fact that was 
found for genus zero in [H] and for genus one in |18j . 

Evidently only for transsonic initial data might the MNLS equation behave unlike either the focusing 
or defocusing cubic nonlinear Schrodinger equation in the semiclassical limit. This observation is the key 
motivation for the study that we begin in this paper. 

To analyze the solution of the MNLS equation ([I]) subject to suitable initial data of the form ^ by 
means of the inverse-scattering transform, one must in general study the small-e asymptotics of two types 
of problems: 

• The direct scattering problem. This amounts to the asymptotic analysis of a singularly perturbed 
linear differential equation (the WKI spectral problem) . The natural tools applicable to this problem 
include the classical WKB method and its generalizations. 

• The inverse scattering problem. This can be formulated as a matrix Riemann-Hilbert problem 
(see ^J3]) with a highly oscillatory or rapidly exponentially growing and decaying jump matrix. The 
natural tools applicable to this problem include the noncommutative steepest descent method of 
Deift and Zhou [7] and its generalizations (in particular the mechanism of the so-called g-function 
first introduced in [DJ). 

Although the asymptotic analysis of the direct scattering problem can be carried out in part using classical 
methods, what really makes the problem difficult is (i) that one requires asymptotics that are uniform with 
respect to the spectral parameter, including near exceptional values where turning points collide, and (ii) 
that for some initial data it is required to approximate some exponentially small quantities that — like the 
semiclassical above-barrier reflection coefficient for the Schrodinger operator — cannot be suitably estimated 
without assumptions of analyticity and the use of turning point analysis in the complex x-plane. For a flavor 
of the type of calculations required to analyze the WKI direct scattering problem for the MNLS equation in 
the semiclassical limit, see [ID] . It should be stressed however, that with very few exceptions in the literatur^] 
the asymptotic information that is available from the direct problem is frequently of insufficient quality 
and/or accuracy to continue with the inverse scattering problem without unjustified formal approximations. 
This is especially true when (as for the WKI problem or the Zakharov-Shabat problem for the focusing 
nonlinear Schrodinger equation) the direct spectral problem cannot be cast as an eigenvalue problem for a 
self-adjoint operator. 

On the other hand, the completely rigorous asymptotic analysis of matrix Riemann-Hilbert problems is 
now a rather well-developed science. In order to bypass the difficulties mentioned above in regard to the 
direct problem and hence enable rigorous analysis of the Cauchy problem for the MNLS equation ([I]), we 
will choose the initial condition functions Aq(x) and Sq(x) to be of the specific form 

(11) ^o^) = ^sech(:c) and So(x) = Sq(0) + Sx + //log(cosh(x)) 

for real parameters v 7^ 0, Sq(0), 5, and ji. The corresponding initial values for the Madelung fields are 

(12) Po(x) = ^ 2 sech 2 (a:) and uq(x) = 5 + /xtanh(a;). 

In fact, without loss of generality, we will assume that v = 1 and So(O) = 0, as this can be accomplished by 
making the substitutions 

(13) <\> = ve iS °(° )/e 4>, e=\u\e, at = \u\- 1 a, t=\v\i 

(x is not scaled) in both ([lj and (111 and dropping the tildes. Unlike the cubic nonlinear Schrodinger 



equation, the MNLS equation is not invariant under Galilean boosts [5] so the parameter S cannot be 
removed even by going to a moving frame of reference. Therefore, both S and [i are essential parameters in 



the family of initial data of the form (111. The advantage of working with the family ( 11 ) is that this choice 



*A notable exception is the analysis of the Korteweg-de Vries equation in the small dispersion limit by Claeys and Grava 
[3] |U [5]. Of course in this case the direct scattering problem corresponds to the well-studied and self-adjoint Schrodinger 
operator, and sufficiently accurate semiclassical approximations of the scattering data for this problem were obtained (only 
fairly recently) by Ramond [20] . 
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renders the direct scattering problem for the MNLS equation equivalent for all e > to a hypergeometric 
equation, and this in turn allows the corresponding scattering data to be obtained without any approximation 
whatsoever [8]. 

For the initial data (12) with v = 1, the quantity Q defined by (10) becomes (for t = 0) 

(14) Q(x) — a 2 sech 2 (a;) + oe5 + a/itanh(x) — 1 = — a 2 T 2 + a/iT + a 2 + aS — 1, T := tanh(x). 

To be in the transsonic case, we need Q(x) to change sign. As a quadratic function of T — tanh(ir) it is 
clear that there are at most two roots. For simplicity we want to arrange that there is exactly one simple 
root of this quadratic in the interval T E (—1, 1), and this will occur if and only if it takes opposite signs for 
T = ±1; hence we assume the condition 

(15) a\n\>\l-a5\. 

For technical reasons that we will explain in §2.2| we also assume the following two conditions: 

(16) (i 2 > 4(1 — aS), 
and 

(17) n > 0. 
If we introduce the combined parameters 



(18) 



A := 



1 — aS 
4a 2 



and B := 



4a' 



then conditions ( 15 )— ( 17 ) imply that (A, B) should occupy the shaded region in the diagram pictured in 
Figure [TJ As an example, some specific parameter values that satisfy these conditions are a = 5 = 1 and 




Figure 1. The shaded region corresponds to the inequalities (15)-(17) 



/i = 2, corresponding to A = and B = 1/2. Under the inequalities (15)-(17), the unique simple root in the 
interval (—1,1) of the quadratic (14) is exactly 



(19) 



T = T c := IB - ^4B 2 -4A+1 



(the other root exceeds T = 1 for (A, B) in the admissible region), and it follows that when t = 0, Q < 
(elliptic case) for x < x c and Q > (hyperbolic case) for x > x c , where x c :— arctanh(T c ). Our main result 
is then the following. 

Theorem 1. Let <f) e (x,t) denote the solution of the Cauchy initial-value problem for the MNLS equation |l]) 
subject to initial data of the form ([2| with Aq(-) and So(-) given by (11), where (without loss of generality) 
v = 1 and iSo(O) = 0, and where a > 0, S, and /i are subject to the inequalities (15)- |l7| ). Then there exists 
a smooth curve x — x c (t), t>0, with x c (0) — x c such that for all t > and all x < x c {t), 



(20) 



Mx,t) = A(x,ty s( - x ^/ e + 0(e), 



0, e > 0, 



where the error term is uniform for (x, t) in compact subsets and where A(x, t) and S(x, t) are smooth, 
real-valued functions independent of e that satisfy A(x,0) — Aq(x) and S(x,0) = So(x). Also, Q < holds 
strictly for all x < x c (t) and t > 0, while Q — > as x — > x c (t), x < x c (t). Finally, whenever x < x c (t) and 
t > 0, the Madelung-type fields 

dS 

(21) p{x,t) := A(x.t) 2 and u(x, t) := — (x, t) 

ox 

exactly satisfy the dispersionless MNLS system ^ . 



The functions A(x,t) and S(x,t) will be specified precisely as part of the proof; see (220). Significantly, 



these functions are obtained without directly analyzing the dispersionless MNLS system, but rather arise 



from the solution of certain algebraic equations (the expressions given in ( 220 ) depend on the solution of ( 79 ) 



and the identities (155) and (160)). Our result is interesting in part because it is a global description of the 
solution <j) e (x,t) valid for all time t > in a semi-infinite spatial domain x < x c (t). Unlike the semiclassical 
limit for the focusing nonlinear Schrodinger equation (for which the dispersionless system analogous to 
Q is elliptic for all initial data) [14j [21] , the solution of the MNLS equation with the initial data under 
consideration does not "break" at some caustic curve beyond which the Madelung ansatz fails and a higher- 
genus multiphase behavior ensues. Indeed, the boundary curve x = x c (t) is merely a sonic line across which 
the asymptotic solution changes from subsonic type (for x < x c (t)) to supersonic type (for x > x c (t)). While 
we do not give the proof in this paper, the sonic transition turns out not to be catastrophic for the Madelung 
fields, which extend smoothly into the supersonic region (although in this region the Madelung theory turns 
out to no longer be globally valid in time, that is, the MNLS solution "breaks" in the hyperbolic region 
for the dispersionless MNLS system). The analysis corresponding to the supersonic region x > x c (t) is the 
subject of our forthcoming paper 9 in this series. 

From the point of view of technique, this paper contains several novelties that we wish to point out. The 
basic approach is to formulate the Riemann-Hilbert problem of inverse-scattering for the MNLS Cauchy 
problem with available explicit formulae for the scattering data (the reflection coefficient is the key quan- 
tity). Being as the formula for the reflection coefficient involves the Euler gamma function with arguments 
proportional to the large parameter e _1 , one can hardly resist applying Stirling-type approximations. While 
much of the analysis corresponds to working with these approximations, we are very careful about not dis- 
carding any error terms until they can be rigorously controlled (which always occurs at the very end of the 
process in the steepest descent method). In particular, note that the function (f> e (x,t) in Theorem [I] is the 
exact solution of the Cauchy initial-value problem for the specified initial data; we do not replace the initial 
condition with a nearby function corresponding to approximate scattering data as is a common approach 
in the subject (see, for example [2J, [T3], and [H]). There are points where the Stirling approximation fails 
(these correspond to the exceptional values of the spectral parameter where turning points collide in WKB 
theory), and we deal with these points with the help of a new kind of "lens-opening" deformation of the 
Riemann-Hilbert problem that completely sidesteps the need for any local analysis near these points. This 
same technique also completely solves certain problems arising from the fact that in cases where one is lucky 
enough that the reflection coefficient admits an analytic continuation away from the real axis, there fre- 
quently exist "phantom poles" that do not correspond to eigenvalues of the spectral problem but nonetheless 
that obstruct the sort of contour deformations that are at the heart of the steepest descent method of Deift 
and Zhou. 

In Sj2]we summarize the key facts about the scattering data for the WKI spectral problem in the case that 



the initial conditions for the MNLS equation are given in the special form (11 ), incidentally correcting a sign 
error in the original reference [5]. Then in ^3] we recall the formulation of the inverse scattering problem for 
the MNLS equation as a Riemann-Hilbert problem of analytic matrix factorization. To prepare the Riemann- 
Hilbert problem for analysis in the semiclassical limit we need to introduce an appropriate "g-function" g(z) 
and a closely related function h(z); these functions are constructed and completely characterized in [j4j We 
should point out that it is the count of the zeros of the auxiliary function Y(z) related to derivatives of g 
and h (see Proposition [T] in SQ that ultimately gives rise to the global-in-time nature of our main result. 
Finally, in f|5]we use the g-function along with the new "lens-opening" method to reduce the Riemann-Hilbert 
problem formulated in Sj3] to a form where it can be treated by means of Neumann series for an associated 
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system of singular integral equations of small-norm type. The estimates resulting from this analysis complete 
the proof of Theorem [I] 

Regarding notation, we use 3?{z} and to denote the real and imaginary parts, respectively, of a 

complex number z. We will typically use the symbol r\ to refer to an arbitrarily small positive quantity. We 
write all matrices in boldface (e.g. M), with the exception of the Pauli spin matrices: 



(22) 



"0 


1" 




'0 -i 




'1 " 


1 





, <J2 ■ = 


i 


, 0-3 := 


-1 



In particular, the expression a 3 appears frequently; this is simply the diagonal matrix 
(23) 

If f{z) denotes a harmonic or analytic scalar or matrix- valued function defined on the complement of an 
oriented arc A in some neighborhood of the complex z-plane, then for z £ A we define f+(z) (respectively 
f-{z)) as the boundary value taken by f(w) asw->zei from the left (respectively right) as A is traversed 
according to its orientation. Throughout our paper, when we write z p for generally complex z and p we 
always mean the principal branch defined by z p := e pi ° s ^ with — n < 3{log(z)} < tt. Finally, we denote 
complex conjugation with an asterisk: z* . 
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a,b E 



2. The Scattering Data and Properties Thereof 

2.1. Exact Formulae. According to the theory of the inverse-scattering transform as described in the 
Appendix of [5] , the key scattering coefficients needed to formulate the inverse-scattering problem are denoted 
S\2{k) (defined for 3{fc 2 } = 0) and S22(k) (defined for 9{fc 2 } > and analytic in the interior of this region). 
For initial data of the special form ([2| with A(-) and S(-) given by (111, the linear scattering problem was 



reduced to a Gauss hypergeometric equation in [8], which allowed the relevant scattering coefficients to be 
calculated explicitly in terms of the Euler gamma function PQ . In the case of v = 1 and Sq — the scattering 
coefficients found in this way are exactljj^] 

(24) SWfc) := - ■ 12 ?' % ] 12 . 2e „ J , 3{fc 2 } = 0, 

and 

V(— _i_ I 1am tv! _i_ H ii±\ 

(25) ^):= r ltnll' }lU + l > ^ 2 }>0- 

Here 

(26) n = fi(Jb) := — ( k 2 - 

ia \ 4 

anc0 

(27) R = R{k) := e l7, ^(-i(16k 2 - ^ 2 )) 1/2 . 

The corresponding reflection coefficient is defined generally in terms of the scattering coefficients as follows: 

(28) r ( k y-=-¥^v k2eR - 

d22\k) 

In Sj3] we will show that we only need to consider k to lie on the boundary of the first quadrant in the 
complex plane (and for some purposes to admit analytic continuation into the open first quadrant), so we 
introduce a new coordinate z as follows: 

(29) z = k 2 and k = e™ l A {-izfl 2 . 



2 Equation [ ]l24[ i corrects a sign error in 8 . In that paper, the sign error originates in equation (146) and subsequently 
propagates into (151), (156), and (158)-(160). 

■^The function R(k) defined in | |27| is actually the Schwarz reflection of the function with the same name used in [8]. 
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Thus, the closed upper half z-plan^] corresponds in a one-to-one fashion to the closed first quadrant in the 
/c-plane. In terms of z, the specific reflection coefficient under consideration is then 

17 W V V ^ ^ 2eW4(_^)i/2 r( i + + r (_ffi _ |)r(-^ + H) 

where fl and i? are written explicitly in terms of z = k 2 . 

A simple exact formula is available for |s(z)| 2 for z £ K as a result of three observations. Firstly, since f2 
is imaginary for z £ M, and since r(w/) = r(iu)*, it follows that 



(31) 

Secondly, the function 



1 \ 2 2e e J 



= f , z e 



TV 4 _i_ n M^wfi 4- 4- iLi 

TV— in _ I^itv— ffi 4- Bl\ 

L \ 2e 2e> L \ 2t^~2e> 



(32) m(z) := 

1 l — 2e 2e^ V 2e 1 2e' 

is meromorphic in a horizontal strip containing the real z axis, and hence so is the function M(z) := 
m(z)m(z*)* . Of course, for z £ E we have M(z) = |m(z)| 2 , and we may calculate M(z) for all real z 
by obtaining a formula valid for z in some interval of the real axis and applying analytic continuation to 
assert that the same formula holds for all z £ K.. For z < y 2 /16, R and SI are both imaginary, so using 
T(w*) = T(w)* we obtain 

f33) M(2) m + " - ire - 1 + mm + g + tni - g - m z< i_ 

[r(-i-i)r(i + i)][r(-^ + |)r(i-i)j ' Me' 

Thirdly, we recall the two reflection identities for the gamma function pQ: 
(34) T(w)T(-w) = r and r(i+w)r(| 



w sin(7rw) cos(7rw) 



Using these in ( 33 ) , we obtain 

M(z) 



(35) 



(zi - iKzi + j) ^(-^ - §22) ^(-^ + &22 

cos(7r(?-|))cos(7r(? + ||)) 
M 2 - j-iR) 2 cosh 2 (%) - cosh 2 ( fg ) 
4e 2 cosh 2 (^) + cosh 2 (§5) - 1 ' 



- 2ia - 

It is clear that the latter is an even function of R and hence defines the analytic continuation of M(z) to the 
whole real z-axis. From these three facts it therefore follows that 



cosh 2 (fi)- cosh 2 (fP) cosh 2 (m- cosh 2 fP), 

36 s z 2 =sgn z „ \ 2i ' „ \_ 2 " ; = 1 „ \ 2e ' „ \_ 2 " jl , z£ 

cosh 2 (^)+cosh 2 (f^)-I cosh 2 (^)+cosh 2 (? R ' 



Since cosh (f^) is even in i?, it follows that the product sgn(z)|s(z)| has a meromorphic continuation from 
K to the whole complex z-plane. 

2.2. Analytic Properties of the Scattering Data. When the scattering coefficients are given by the 



specific formulae p4[)-(25), it is obvious that S\2 and S22 are even functions of R and hence they have no 
branch point at z — fi 2 /I6, the square-root branch point of R considered as a function of z — k 2 . This makes 
kSi2(k) and S22 (k) both meromorphic functions of z = k 2 . From this it follows that s(z) has a meromorphic 
continuation from the real z-axis into the full upper half z-plane. The point z = is a branch point of s(z), 
however, as is evidenced by the fact that the product s(z)(— iz) -1 / 2 is analytic and non- vanishing near z = 
(the apparent singularity at z = from the factor (— iz)~ x / 2 is cancelled by the simple zero of T(— ^ + ||) _1 ). 
This shows that s(z) does not have any single- valued continuation into the lower half z-plane due to the 
mismatch of boundary values taken on the branch cut of (— iz) 1 / 2 along the negative imaginary axis. In this 
paper, we will refer to both the meromorphic continuation of s(z) into the upper half z-plane as well as its 



The coordinate 2 denned in | |29[ l is the negative of the variable 2 used in [8] . 
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boundary value taken on the real z-axis by the same notation: s(z). Now we consider the analytic nature of 
s(z) for > 0. 

In general, any zeros of the scattering coefficient S l 22(fc) for 3{fc 2 } > (equivalently for > 0) 

represent discrete spectrum (eigenvalues) of the direct scattering problem and ultimately poles in the matrix 
unknown of the Riemann-Hilbert problem of inverse-scattering theory. For technical reasons we would like 



to avoid having to include such poles. This motivates the condition (16) on the parameters a, S, and fi; 



indeed in the case that 822(h) is given by the formula (25), it was shown in [5] that under the condition (16) 
there are no zeros in the closed upper half z-plane for any e > 0. In this case, 6*22 is an analytic function of 
z for > that extends continuously to the real z-axis and that for each e > is bounded away from 

zero. 

The scattering coefficient 6*12 generally has no analytic continuation from the axes 3{fc 2 } = (or the real 
axis Q{z} = 0). However, in the case of the special formula (24), it is clear that S12 can be continued into 
the upper half z-plane as a meromorphic function. The factors r(| — ^ + ^) and r(— ^ + |^) are analytic 
and nonvanishing for > 0, but the factor r(| — ^ — — ) contributes to S±2 an infinite array of simple 

poles with small spacing ae/2 along the vertical ray 3?{z} = zp and 3{z} > where 



(37) 



aS + a/j,). 



These poles are singularities of s(z) that unlike zeros of 622 do not correspond to eigenvalues of the direct 
scattering problem; in [8 they are called phantom poles. The presence of such poles in the reflection 
coefficient of an analogous inverse-scattering problem for the focusing nonlinear Schrodinger equation has 
been previously dealt with [211 by means of the installation of an implicit local parametrix; by contrast 
the deformations of the Riemann-Hilbert problem that we will carry out in <|5] are specially designed to 
render this unfortunate feature of the reflection coefficient completely harmless without the use of any local 
parametrix at all. The remaining factor r(— ^ — ^) can only contribute zeros to Si2- For technical reasons 
we wish to prevent these zeros from lying in the upper half z-plane, which is the purpose of the condition 



(17) on /i. Were this condition not satisfied, the factor in question would contribute to S12 infinitely many 



simple zeros lying along (and densely filling out as e 4- 0) a parabolic curve in the upper half-plane, and their 



presence would cause difficulties with the use of Stirling's formula for asymptotic analysis (see [2.3). 

Under assumptions (16) and (17), s(z) is analytic and non-vanishing in the open upper half-plane with 
the exception of the "phantom poles" along the line 3?{z} = zp. Note also that the condition (16) implies 
that the zero z = ^t 2 /16 of R is greater than the zero z = (1 — aS)/A of f2, and that conditions ( 15 ) and ( 17) 
taken together imply that zp > 0. 

2.3. Asymptotic Properties of the Scattering Data. 



2.3.1. Asymptotic behavior of \s(z)\ 2 . The three conditions ( 15 )— ( 17 1 also determine the qualitative asymp- 
totic behavior of |s(z)| 2 given by (36) for zGMin the limit e I as we will now show. We have O G iR for all 
zel, and jjL is a real constant. The relative size of the positive real (quadratic and constant, respectively) 
quantities 



(38) 



q(z) := Hft) 2 



EL 
4 



J-^(4z-(l-^)) 2 and c(z):=(£ 

will therefore play a role in determining the asymptotic behavior of |s(z)| 2 for real z in the limit e l 0. If 
z > /i 2 /16, then R £ R and these two functions determine the asymptotic behavior completely; on the other 
hand, if z < /^ 2 /16, then R £ iR, and we also need to take into account the positive real (linear) quantity 



(39) 



l(z) :-- 



.R 



EL 
4 



4s, 



z < 



EL 

16' 



Extending the definition of l(z) by setting l(z) := for z > /i 2 /16, we observe that the asymptotic behavior 
of |s(z)| 2 is determined by the maximum of q(z), l(z), and c(z): 

• For those z for which c(z) = m&x{q(z),l(z),c(z)}, |s(z)| 2 is exponentially large in the limit e — >• 
(exponential growth of s(z)). 

• For those z for which l(z) — max{<7(z), l(z), c(z)}, |s(z)| 2 — 1 is exponentially small in the limit e — > 
(pure oscillation of s(z)). 



• For those z for which q(z) = max{g(z), l(z), c(z)}, |s(z)| 2 is exponentially small in the limit e — > 
(exponential decay of s(z)). 

It is easy to see that for z < 0, l(z) > c(z), and that for z > 0, l(z) < c(z), while Z(0) = c(0). It follows that 
there are only three possible scenarios for how the real z-axis may be partitioned into intervals of growth, 
pure oscillation, and decay of s(z), as illustrated in Figure [2] Noting that the points z = (1 — a6 ± a|/i|)/4 




Figure 2. The three configurations of the reflection coefficient. Left: 1 — a5 + a\/i\ < 0. 
Center: a\n\ > |1 — a6\. Right: 1 — aS — a\fi\ > 0. The subdominant arcs of the graphs of 
q(z), l(z), and c(z) are shown with dashed curves. 



are the solutions of the equation q(z) = c(z), these scenarios are: 

• If 1 — a5+a|^i| < 0, then there are two real and negative roots of q(z) = l{z), and for z between these 
two roots s(z) is purely oscillatory, while for z outside of this interval s(z) is exponentially small. See 
Figure [2j left-hand panel. In such cases one expects the asymptotic analysis of the Riemann-Hilbert 
problem of inverse scattering to resemble that which has been carried out for the Korteweg-de Vries 
equation in the small dispersion limit (see [5] as well as [3]-[5]). In such problems all of the necessary 
deformations of the Riemann-Hilbert problem that are required for the steepest descent technique are 
purely local to the real axis in the z-plane, and hence analyticity of the scattering data does not play 
any central role in the theory. This condition really corresponds to the case of globally supersonic 
initial data for the MNLS equation, for which there is a clear analogy with the inverse-scattering 
theory of the defocusing nonlinear Schrodinger equation in the semiclassical limit as shown in |10) . 

• If 1 — aS — a\n\ > 0, then s(z) is exponentially large for (1 — a>8 — a\/i\)/4 < z < (1 — aS + a\/i\)/A 
(where c(z) is dominant) and otherwise is exponentially small. See Figure[2j right-hand panel. In such 
cases one expects that to handle the exponentially large reflection coefficient it will be necessary to 
introduce the meromorphic continuation of s(z) from its interval of growth into the complex z-plane. 
This type of analysis was carried out for the semiclassical limit of the focusing nonlinear Schrodinger 
equation in |21j . In fact, this condition corresponds to the case of globally subsonic initial data 
for the MNLS equation, and the analogy with the focusing nonlinear Schrodinger equation in the 
semiclassical limit is also considered in |10) . 

• If a\fi\ > \1 — a5\, (that is, 
by z = zl := z~, where 



(40) 



if (15) holds) then there is a unique negative root of q(z) = l(z), given 



(1 - aS) 



a 2 ± 



a 2 fi 2 



a 2 (l — a5) 



and s(z) is exponentially small for z < z^ and for z > 
zl < z < 0, and is exponentially large for < z < (1 - 



(1 — aS + a|/i|)/4, is purely oscillatory for 
- aS + a|/x|)/4. See Figure [2j center panel. 
Note that under the assumption ( 17), the point of transition from exponential growth to exponential 
decay is exactly z = zp. 

which was originally imposed to ensure that the initial condi- 



These results show that the assumption ( 15 



tions were of transsonic type, is also exactly what is required to ensure that in the spectral transform domain 
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there exist both a negative interval of pure oscillation and also an abutting positive interval of exponential 
growth of the reflection coefficient s(z). 

To handle the Riemann-Hilbert problem of inverse scattering in the transsonic case under consideration will 
therefore require in the same problem a combination of techniques from both the category of "modulationally 
stable" semiclassical limits (e.g. Korteweg-de Vries, defocusing nonlinear Schrodinger) and also the category 
of "modulationally unstable" semiclassical limits (e.g. focusing nonlinear Schrodinger). 



Assuming the three conditions ( 15 1 — ( 17 ) , note that as shown in the diagram in the central panel of 



Figure [2] there are five distinguished points on the real z-axis, in order from left-to- right: 

• z = zl, the point of transition from exponential decay to oscillatory behavior of s(z). This is the 
negative root of the quadratic equation q{z) = l(z). 

• z — •= j(l — oc5 — a/i), a point at which two subdominant exponentials in |s(z)| 2 exchange roles. 
This point is the negative root of the quadratic equation q(z) = c(z) and it also has significance with 



respect to the turning point curve to be explained in S4.1 

• z = 0, the point of transition from oscillatory behavior to exponential growth of s{z). This is the 
unique root of the linear equation l{z) — c(z). 

• z = z + , another point at which two subdominant exponentials in |s(z)| 2 exchange roles. This point 
is the positive root of the quadratic equation q(z) = l(z) and is given explicitly by (40). 

• z = zp :— |(1 — ad + a/i), the point of transition from exponential growth to exponential decay of 
s(z). This is the positive root of the quadratic equation q(z) = c(z). 

That z ± are real and distinct follows from condition (16 1 along with a > 0, and that they have opposite 



signs is exactly equivalent to the condition (15 1. Note also that conditions (151 and (17 1 together imply that 
z + < zp. For later convenience, let us define the linear exponent /(z) by 



(41) 



/(*) 



2tt 



(z - z P ). 



The assumptions (151 - (17) imply that simple bounds for |s(z)| 2 then follow from the exact formula (36) 
and the central graph in Figure [2] Indeed, since cosh 2 (7ri?/2ie) > cosh 2 (7r/i/2e) and cosh 2 (7rfi/«e) — 1 > 
for z < 0, we easily obtain the inequality: 



(42) 



|s(z)| 2 < 1, z<0. 



On the other hand, for z > we have instead that 
\s(z)\ 2 < 



2 ^ cosh 2 (^) 



cosli (^) 



(43) 



_ e 4 7r (a M /4-|^-(l-Q5)/4|)/(Q£) 



< 4 e 4^(ap/4-|2-(l-a<5)/4|)/( Qe ) 

= 4e- 2 ^>/ e , z>0. 



1 + e-^v/t 

I + e -47T\z-(l-a8)/4\/(ae) 



Note that by the Mean Value Theorem applied to ( 36 1 , the estimates ( 42 ) and ( 43 ) can be replaced by an 



improved estimate valid in a neighborhood of the origin as follows: for each sufficiently small 77 > there 
exist constants C > and K > independent of e such that 



(44) 



C\i 



\z\ < rj. 



Here we have used the fact that zp > to ensure that K > 0. 

Similar reasoning produces the following asymptotic formulae involving |s(z)| 2 . 



(45) 1-Kz)| 5 



-T(»/e 



(1 + exponentially small in e) 
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uniformly for z € [zi 



Firstly, we have 

V] U [zoo + V, 0] 



for all i] > 0, where t(z) > is defined by 

r ZL (z) := -mR - 2tt|£1| = ttvV - 16z + 



(46) 



r(z) 



4tt 

Q 



1-aS 



Z L < z < z a 



(,Tq(z) := —ITTR — TTfl = -K\J fl 2 — 16z — TTfl, z^^zkO. 

The function t(z) extends (by zero) to a continuous function for z £ K. Secondly, we have 
(47) l s ( z )| 2 = e~ 2 -"^/ e (l + exponentially small in e), uniformly for z > z + + 77, \/rj > 0. 



2.3.2. Stirling asymptotics for s(z) . Under conditions ( 15 1— ( 17 1, Stirling's formula [T] yields accurate asymp- 
totics for s(z) that are uniformly accurate to a relative error of order O(e) in the closed set S consisting of 
the closed upper half-plane with the following sets omitted: 

• The vertical strip — Zp\ < r\ with > 0. This strip contains the accumulation locus of the 
phantom poles. The Stirling approximation of the ratio r(| — ^ — 7)/!" (^ + 7^ + 7) fails (only) in 
this strip as e 1 0. 

• The open upper half-disk of radius 77 centered at z = 0. The point z — is the unique root of R—ifx 
for Q{z} > when fi > (and R + i^j, is nonvanishing for > 0). The Stirling approximation of 
the product 



Y t ) fails (only) in this half-disk as e l 0. 



• The open upper half-disks of radius 77 centered at the real points z~ < < z + defined by (40 1. 
These are the roots of the quadratic equation R 2 — (2f2) 2 . The Stirling approximation of the 
product r(| + 7 - 5j)r(^ + 7 + J|) fails (only) in the union of these half-disks as e | 0. 

Here 77 is any arbitrarily small fixed value. The omitted vertical strip containing the phantom poles divides 
S into two disjoint subregions: Sl to the left of the strip and Sr to the right. In each of these two subregions 
Stirling's formula will supply a different analytic approximation to the ratio s(z). The resulting formulae 
are as follows. Firstly define (here log(-) denotes the principal branch with |3{log(-)}| < it) 



9{z) 



4 (1 
-fi log(2) + - (z - zp) log - (zp - z j 



(48) 



iR 



2 

CI - 



log 
R 



—ift — R 



fi 



iR 



log n - 



R 



2 

- i (n 



log 

R 

^2 



-ifi 



R 



2 

log (fi 



> 0. 



Under the conditions (15|-(17|, <&(z) is an analytic function of z for 3{z} > 0. Now, define: 

-,s(z)e~ ltE,(z)/e , for > 0. 



(49) 



E(z) :- 



The function E(z) is meromorphic for > and its poles are confined to the line 3?{z} = zp. In 

particular it is analytic in the interior of S = <Sl U Sp. Stirling's formula implies that 



(50) 



E(z) 



l + 0(e), zeS L , 
e -2/W/« (i + 0(e)) , z e S R , 



where f(z) is defined by (41), with both estimates holding uniformly in the indicated region (including for 
the boundary value E + (z) taken on the real axis). Next, define 



(51) 



E{z) := (l + e 2/ W/ e )£(z), > 0. 



The explicit prefactor exactly cancels the poles in E{z) for > 0, and hence E(z) is an analytic function 

in the whole upper half complex plane. Letting S denote the closed upper half-plane with the open upper 
half-disks of radius r\ centered at z = zp, 2 = 0, z = z + , and z — zp omitted, we can obtain from Stirling's 



formula (by using first the reflection identities (34 1 for T(-)) that 



(52) 



E(z) = 1 



O(e), zeS. 
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Again, the estimate holds uniformly (including for the boundary value E + (z)). It will be important later to 
record the imaginary part of the boundary values taken by <£> on the real axis: 



(53) 



(R(z) - 2Sl(z)) , 



3{$+(z)} = ^ 



ITT 

■— {R(z) - in) , 



2tt 



a 



\z- Zp\ 



-1/001. 



z < Z L 



Z\ J = z < z < 



< z < z H 



Z > Z n 



It therefore follows that &(z) extends through the interval zl < z < to an analytic function defined 
also in the open lower half-plane, and henceforth we consider $>(z) to be analytic in the slit domain z g 
C\((— oo, zl]U[0, +oo)) and to satisfy the Schwarz symmetry condition $(z*) = <I>(z)*. By direct calculation, 
one sees from (|48l) that 



(54) 



2vri 



-i/2 _ a. 



0{z 



-3/2 n 



oo, 3{z} > 0. 



An important observation is that while here the formulae (46) and (48 1 for t(z) and $(0) respectively 
come directly from analysis of the exact formula ( 30 1 for the reflection coefficient valid in the case of the 



particular initial conditions under consideration, these formulae agree exactly with predictions based on 
WKB theory formally valid for more general initial data. This formal semiclassical spectral analysis of the 
direct scattering problem can be found in all details in §2 of |10j . The more general formulae corresponding 



to (46 ) and (48 1 arise from WKB theory as follows. The turning point curve of WKB theory is, by definition, 



the locus of real x and complex z satisfying the characteristic equation 

(55) x(x; z) := {Az - 1 + au (x)) 2 + 16a 2 zp Q (x) = 0. 

Let z € K be a real value for which the function x( x i z ) nas re& l roots x, and let x_(z) denote the most 
negative of these roots. In particular, x(x; z) > for all x < X-(z). Define signs a± by 

(56) a± := sgn(4z — 1 + au±), 
and then set 

(57) w(x;z) := cr_ \/ x(x; z), x < X-(z). 

For each appropriate value of z € K, this quantity has the fixed sign cr_ throughout its domain of definition, 
and uj(x; z) — > ui-(z) asn-> —00, where 

(58) oo- (z) := Az - 1 + «m_ . 
Then, for all real z for which x(x; z) has real roots, we let 

(59) 

and then 



7(2; z) := y/-x{x;z) > 0, x( x \ z ) < °: 



(60) 



-{z) = 



1 



is a generalized formula that agrees exactly with (461 when Aq(x) and So(x) are the particular functions 



(11) considered in this paper. Next, we have 

I /•*-(*) 1 
(61) $(z) = - / [u(y;z)-u-(z)] + -w_(z)x_(z) + 5_ 



This formula is a generalization of ( 48 ) in the sense that the latter is the analytic continuation to the upper 



half z-plane of (61 1 in the special case of the functions Aq(x) and Sq(x) given by ( 11 ). For these functions 



it turns out that there are generically either two distinct roots (turning points) X-(z) < x + {z) of x{ x \ z ) 
or none for z € K. In such cases it can be shown that if zq is a real value at which x±(z) coalesce in a 
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square-root sense, then the functions $(z) and r(z) can be analytically continued about such a point, and 
the following monodromy relations hold: 

(62) $(z + - z )e ±2m ) = $(z) ± ia-r{z), 



(63) 



r(z + (z- z )e ±2m ) = r(z) 



We will make use of the concrete formulae ( 46 1 and ( 48 1 and also the equivalent WKB formulae ( 60 ) and 



(61 1 at various points in our analysis. 

3. RlEMANN-HlLBERT PROBLEM OF INVERSE SCATTERING 

3.1. Basic problem. We take as a starting point the formulation of the inverse-scattering problem as 
a matrix Riemann-Hilbert problem for an unknown *M.(k:x,t) as described in the Appendix of 8J. This 
matrix has jump discontinuities across the real and imaginary axes in the complex fc-plane and satisfies the 
holomorphic involution M(— fc; x, t) = i CT3 M(fc; x, t)i^ rT3 and the antiholomorphic involution M(— k*; x, t)* = 
criM(fc; x, t)<T\. If for 3{fc 2 } = 0, M±(fc; x,t) denotes the boundary value taken by M from a quadrant 
where ±Q{fc 2 } < 0, then 

i6(k 2 ;x,t)a 3 /t I" 1 ± \ r ( k )\ 2 r ( k ) 
±r(k)* 1 



(64) 
where 
(65) 



M + (k;x,t) = M_(k;x,t)e 



-i6{k 2 ;x,t)a 3 /e _j_^,2 ^ g 



6{z) = 9(z;x,t) := 





4 




h: 


X 2 






or 





and where r(fc) 

S 12 Q 
(66) 



-r(—k) is the reflection coefficient defined generally in terms of the scattering coefficients 
Si2(k) and 5*22 (k) by (28). The matrix M(fe; x, t) is normalized to the identity at k = 0: 



lim M(k;x,t) = I 

k—tO 



regardless of the quadrant from which the limit is taken, and is required to have a well-defined limiting 
value (a matrix- valued function of (a;, t)) as k — > oo. In the case under study in this paper that there are 
no eigenvalues of the scattering problem, the matrix M(fc;x,i) is analytic for 3{/c 2 } ^ and the above 
properties uniquely determine M(fc;x, t) given r(-). The corresponding solution of the MNLS equation is 
then obtained from 'M.(k;x,t) by the formula 

2k M 12 (k;x,t) 



(67) 



4> e (x, t) := lim 



An-oo a M22(k;x,t) 

(Note that M(fc;x, £) is asymptotically diagonal as fc — > oo. See [51 page 990].) 

3.2. Symmetry reduction. It will be convenient to take advantage of the fc H> — fc symmetry to define a 
new unknown by setting N(z;x,t) := M(z(— z) 1 / 2 ; x, t) for ^ 0. Since (— z) 1 / 2 denotes the principal 

branch of the square root, the matrix M(fc;a;,t) as defined in the first (second) quadrant of the fc-plane is 
equivalent to the matrix N(z;x,t) in the upper (lower) half z-plane. For z£l, let 'N±(z;x,t) denote the 
boundary value of N taken from C±. Then, if z > 0, since N±(z;x,t) — Mz F (±z 1 / 2 ; x, t) and according to 
the holomorphic symmetry of M(fc; x, t) we have M + (— z 1 / 2 ; x, t) = i ' 3 M + (z 1 / 2 ; x, t)i~ a3 , 



(68) N+(z;x,t) =i a3 N_(z;x,t)r' j3 e ie(z '' x ' t) ' j3/e 



-riz 1 ' 2 ) 
l + |r(^/2)|2 



1 

. r ( 2 l/2 

On the other hand, if z < 0, then ~N±(z;x,t) = Mzp(j(— z) 1 / 2 ; x, t), so 

1 -r{i{-zf/ 2 ) 



-i0(z;x,t)(T3 / e 



-i9(z;x,t)(T3 1 e 



z > 0. 



z < 0. 



(69) N + (z;x,£) = N_(z; a; ,i)e^^ t )-/ £ y^^y x _ | r ( i( _ z )i/2)|3| 

The function e l7r / 4 (— iz) 1 / 2 (principal branch) is analytic in the upper half z-plane and agrees with z 1 / 2 for 
z > and with i(— z) 1 / 2 for z < 0. Note that both r(z 1 ^ 2 ) for z > and r(i(— z) 1 / 2 ) for z < correspond 
to the same function s(z) defined for all z G R by (30 1. We frequently omit the explicit dependence on 
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the parameters (x,t) (along with e, a, S, and /i) and simply write N(z) = N(z;x,i). The matrix N(z) 
necessarily satisfies the Schwarz symmetry property: 



(70) 

as well as the normalization condition 
(71) 



N(z*) = o-xN^Vi 



lim N(z) = I, 

z-yQ 



and we require N(z) to have a well-defined (necessarily diagonal) limiting value as z — > oo. 

4. The Functions g(z) and /i(z) 

In this section, we construct two related analytic functions of z that will depend parametrically on 
(x,t) € M 2 . For notational convenience, we will frequently not write the parameters explicitly as arguments 
of these functions, but the reader should be aware that the dependence is there. We wish to stress from the 
outset that all objects considered in this section are independent of the basic asymptotic parameter e > 0. 

4.1. The WKB turning point curve. In the case of the special initial data under consideration the 



turning point curve defined in general by (551 is given by the equation 

(72) (Z-A + BT) 2 + Z(1-T 2 ) = 0, z = a 2 Z, T := tanh(x), 



where the parameters A and B are defined by (18 1. This is quadratic in Z for fixed T and also is quadratic 
in T for fixed Z. In solving for Z, one easily finds that the discriminant is -(1 - T 2 )(T 2 - 4BT + AA - 1), 
whose second factor is proportional by a negative constant to the quantity Q at t = as defined by (14 1. 
Therefore under the conditions (15)-(17), the only (simple) root in the range \T\ < 1 is T = T c defined 
by (19). (The other root of the quadratic factor exceeds T = 1 for (A,B) in the admissible region shown 
in Figure [I]) Moreover, the discriminant is positive (only) for T G (T c , 1). The point x c — arctanh(T c ), 
corresponding to z = z c := (a(S + fiT c ) — l)/4 = (auo(x c ) — l)/4, is the only point on the turning point 
curve in the (x, z)-plane where there is a vertical tangent. For x > x c the turning point curve has two real 
branches that we denote as z — a(x) and z = b(x) with a(x) < b(x), with equality holding only at x = x c 
and in the limit x — > +oo as illustrated in Figure [3] 



Horizontal tangents to the turning point curve correspond to simultaneous solutions of ( 72 ) and the 
identity 



(73) 



B(Z- A + BT) - ZT = Q 



obtained from (72 1 by differentiating implicitly with respect to T and setting dZ/dT = 0. Using this latter 



(T-B) 
A/B or T = T± (A, B) where 

T ± {A,B) :=±- 



2 [{B 2 - A)T 2 + BT - B 2 ] =0 



relation to eliminate Z from ( 72 ) we obtain an equation for T- values of horizontal tangents in the form 

BT-A 

(74) 

Therefore, either T 
(75) 



2B 



V4B 2 -4A + 1±1 

Since B > and B 2 > A both hold throughout the domain of values of (A, B) consistent with conditions 



([T5])-(|T7]) (see Figure [TJ, it is clear that T + (A,B) > while T_(A,B) < 0. Moreover, T + {A,B) = 1 implies 
that B — A while T^(A,B) = —1 implies that B = —A, both conditions that only occur on the boundary 
of the admissible region for (A,B). Also, from the asymptotic expansions of T±(A,B) in the limit B j" oo 
with AeE fixed one concludes that < T + (A, B) < 1 and T-(A, B) < -1 both hold in the interior of the 
admissible region, and therefore of T± (A, B) only T+(A, B) corresponds to a value of x E R. Since B > \A\ 
holds in the interior of the admissible region, T = A/B also corresponds to a real value of x. Similar analysis 
proves that both T = A/B and T = T+(A, B) exceed T c for (^4, B) in the admissible region, so only the real 
part of the turning point curve for T > T c has any horizontal tangents (defined by dZ/dT = 0), and there 
are exactly two of them. Also, it is easy to check that throughout the admissible region of the (A, £?)-plane, 
the horizontal tangent at T — A/B occurs on the more positive of the two branches of Z(T) (with value 




Figure 3. The real part of the turning point curve for (A,B) in the admissible region 
described by conditions ( 15 )— ( 17 ). The red part of the curve is the graph of a function 
x = X-(z) and the blue part of the curve is the graph of a function x — x + {z) both over the 
interval z^ < z < 0. For such z, x±(z) are real turning points of WKB theory. The upper 
curve is the graph of a function z = b(x) and the lower curve is the graph of a function 
z = a(x), both defined for x > x c . The singularity of x+{z) occurs when Z = A — B or 
z = : = (1 — ctS — a/x)/4, the break point of the equation (46 1 for t(z) in the interval 
z L < z < 0. 



(76) 



, i A 

x\ :— arctanh ( — 
B 



and x-x := axctanh(T+(A, B)). 



Note that under the conditions ( 15 1— ( 17 ) , we have the inequalities Xi > X\ > x c . It is also possible to show 
that these conditions imply that z c > z x := (1 — aS — afi)/4 (z^ as originally defined in {2.3 is also the 
point at which the branches of the turning point curve coalesce in the limit x — > +oo, or equivalently, the 
singularity of x+(z)), with equality occurring only along the part of the boundary of the admissible region 
of the parameter space illustrated in Figure [l] with A > 0. This in turn implies that there exists exactly one 
value of x > Xi that we denote by x^, at which point b(x^) — z c . For x c < x < X3 we have a(x) < z c < b(x) 
and for i > i 3 we have a(x) < b(x) < z c . It also follows from the conditions (15)-(17) that in fact £3 > x 2 - 
These relationships are illustrated in Figure [3J 

For x < x c , the turning point curve is complex, with two distinct complex-conjugate branches: z = $(x) 
with 3(3(2;)} > and z — $(x)* . The two branches coincide (with a real value) only in the limits i|i c 
(with value z = z c ) and x \, —00 (with value Z = A + B>0otz — zp). Moreover, it is easy to prove that 
for (A, B) in the admissible region, $t{$(x)} is monotone decreasing in x for — 00 < x < x c , and 



(77) 



sup $tti{x)} = 

x<x c 



lim dl{%(x)} = zp. 



In particular, ?R.{%(x)} < zp for all finite x < x c . Note also that for x < 



(78) 



= (au (x) - 1 + 2a 2 p (x)) . 



1.1 



4.2. Basic construction of g(z) and h(z). The auxiliary function Y(z). Let a complex number q with 
^{q} > be given, along with a simple oriented arc B from z — q to z = zp in the open upper half-plane. 
Let S(z) = S(z; q, q* , B) be the function uniquely defined by the following properties: 

• S{z) is defined and analytic for z £ C \ (B U B*). 
. S(z) 2 = (z- q)(z-q*). 

• S(z) = z + 0(1) as z ->• oo. 

At times we will want to think of q and q* as being independent complex variables, but when they are linked 
by complex conjugation, S(z) is obviously a Schwarz-symmetric function: S(z*) = S(z)*. Note that S(z) 
changes sign across its branch cut B U B* . 

Let us assume (this will be completely clarified in E 4.3 ) that given (x,t) £ R 2 , the complex number q 
satisfies the equations 

2ir 

;x,t) := — (2(q + q*)t + ax-t)+I {q,q*) = 



(79) 

where 
(80) 



M (q,q* 
Mi(q,q* 



or 

IT 



X,t) := — ((3q 2 + 2\q\ 2 + 3q* 2 )t + (q + q*)(ax - t)) + h(q, q*) = 0, 



$'(s) S p ds 



p = 0,l. 



As usual, the notation S+(z) denotes the boundary value taken from the left side of B as the arc is traversed 
according to its orientation from q to zp. Now, let p(z) be the function analytic for z G C \ (B U B*) given 
by the Cauchy-type integral formula 



(81) 



p(z) 



S(z) 
2m 



$'(s) ds 
S+{s){s - z) 



?'(z) 



8tS(z) 



zeC\(B\JB* 



where the notation (*)* indicates the Schwarz reflection w(z*)* of the function w(z) immediately preceding 
the minus sign, and where 8(z) = 9(z; x, t) is defined by (|65[). The conditions (79) obviously imply that 



(82) p(z) 
and hence the contour integral 

(83) g(z) :-- 



z —> oo, 



p(s)ds, zeC\(BUB*) 

is independent of path (as long as the path avoids the branch cut B U B* of p) and defines a function analytic 
in the same domain as p, with a well-defined limiting value 5(00). In terms of g(z) we define a related 
function h as follows: 



(84) 



h(z) := 6(z) + hb{z) 



- g(z), zeC\ ((-00, z L ] U [0, +00) U6U B*). 



Because 9(z) is entire, h(z) is analytic exactly in the intersection of the domains of analyticity of $(z) and 
g{z) as explicitly indicated. Both g and h are Schwarz-symmetric functions: g(z*) — g{z)* and h(z*) — h(z)* . 

Let the function Y(z) = Y(z;q,q*) be defined for > and z near B by the following Cauchy-type 

integral: 

f <f>'(s)ds 
f L S(s)(s - z) 



(85) 



Y(z;q,q*) 



1 

2tti 



a 



where L is a clockwise (negatively) oriented loop in the upper half s-plane beginning and ending at zp and 
encircling both s = z and the arc B exactly once. Obviously Y is an analytic function of z in the interior 
of the loop L. Exploiting more detailed information about $ (as is available from the specific formula (48) 
but that may or may not be available for more general initial data when one uses instead the WKB formula 
(61 )) allows us define a global analytic continuation of Y(z; q, q*) as we will now show. Indeed, using the fact 
that $'(z) is analytic in the open upper half-plane and satisfies <f>'(z*)* — &'(z), along with the asymptotic 

16 



estimate <&'(z) = 0(log |z|) as z — > oo allows us to deform the contour L to a contour along the real axis; 
therefore 



(86) 



Y(z;q,q*) 



16t 

^2 



1 

2?ri 



defines the analytic continuation of Y to the domain z G C\ 



16t 

O 2 



5( S )( S -z) 



c/.s 



Now, we recall ( 53 ) along with the definitions 



y/(s - 5?{g}) 2 + 3{g} 2 for s < z P and 



([26]) of fi(z) and (|27) of and the fact that for s G E, S(s) : 
S(s) = ^(s - 3?M) 2 + 3{ 9 } 2 for s > zp, and we therefore find that when $ is given by (J48J) , 



let 

7a 



(87) r(z;g,<f): 

2 r ZL 



c/.s 



%}) 2 + 3{<?} 2 (s- z) 

f/.S 



(is 



" 16V( S - ^M) 2 + »{q} 2 (s - z) 



vV - 16V(« - 3*{9}) 2 + »M 2 (i 
2 /-+ 00 ds 
a 



+ V(*-»{«}) 2 +S>{9} 2 («-«) 



This completes the extension of Y{z; q, q*) as an analytic function from z near B to the maximal slit domain 

C+UC_U(z L ,0). 

The function Y(z; q, q*) has many purposes in our analysis. For example, it is easy to check with the use 
of elementary contour deformations that the following identities hold: 



g' + (z) - g'_{z) = S + (z)Y(z;q,q*), z e B, 
and 

(89) 2h'(z) = -S(z)Y(z;q,q*), z G C \ ((-oo, z L ] U [0, +oo) U B U B*). 

The key property of the analytic function Y(z;q,q*) that we will require frequently is the following. 

Proposition 1. For t > there exists exactly one simple zero of Y(z; q, q*) in its domain of analyticity 
(and hence by Schwarz symmetry necessarily located in the real interval zl < z < 0). For t < there exist 
exactly three zeros (counting with multiplicities) ofY(z;q,q*) in its domain of analyticity, and in this case 
given any M > 0, for —t sufficiently small exactly one simple zero lies in the real interval zl < z < while 
there is an additional simple zero in each of the upper and lower half-planes with \z\ > M . (In the boundary 
case of t = the non-real zeros are at z = oo.) 

Proof. The domain of analyticity of Y(z; q, q*) is the slit domain z G C \ ((— oo, zi] U [0, +oo)), and we note 
that Y as given by (87) satisfies the Schwarz symmetry condition Y(z*;q,q*)* = Y(z;q,q*). From these 
facts and the Plemclj formula we find that upon taking a boundary value as z approaches the real axis from 
the upper half-plane, 



(90) 



Ana 



16z 



ct\[p? 



0, 



16z^(z-5R{ g }) 2 +9{ 9 } 2 



%{Y + (z;q,q*)}= { 



-47T 



V/x 2 - 16z y/(z - 3?M) 2 + 3{?} 2 ' 



2tt 



z < Z L 



z L < z < 



< z < z H 



z > z n 



R{«}) 2 + 9{ ? } a ' 

Also, from d87| it is clear that the complex boundary value Y + (z;q, q*) is actually analytic except at the 
points z G {zl,0,z + }, where it becomes infinite in magnitude with logarithmic singularities. We may 
combine this information to show that, when the parameters correspond to the admissible region shown in 
Figure [T] the following statements are valid. 
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(1) There is a positive constant c > such that as z increases from z = to z = z + , ^s{Y + (z; q, q*)} < 
-c < 0, while $t{Y + (z; q, q*)} -> -oo as z J, and 3?{F + (z; g, q*)} -> +oo as z t z+. 

(2) As z increases from z = z + to z = +oo, 3{Y+(z; g, g*)} > with equality only in the limit z j" +oo, 
while 5i{F + (z; g, g*)} +oo as z | z+ and 5i{F + (z; g, g*)} 16i/« 2 as z f +oo. 

(3) As z increases from z = — oo to z = zl, g, g*)} < with equality only in the limit z J, — oo, 
while 3?{F+(z; g, g*)} 16i/a 2 as z | -oo and 5i{F + (z; g, g*)} +oo as z | z L . 

In particular, as z increases from z = to z = +oo and then from z = — oo to z = zl, the complex 
boundary value K+(z;g, q*) traces out a curve in the K-plane that crosses the real axis only at one finite 
point, Y = 16t/a 2 , corresponding to z — ±oo. The three intervals in the z-plane and their images under Y + 
in the K-plane are shown in red, with numbers corresponding to the above enumeration, in Figure [4] The 




FIGURE 4. The domain of analyticity of Y(z;q, q*) is C \ ((-oo,z L ] U [0,+oo)), and a 
qualitative sketch of the image of the boundary of this slit domain in the K-plane. 

desired result then follows from an elementary application of the Argument Principle using the information 
enumerated above regarding the boundary value Y" + (z;g, q*) on the cut (and corresponding information 
regarding the boundary value Y—(z; g, q*) follows by Schwarz symmetry). It is obvious that as the boundary 
of the slit domain is traversed once in the positive sense (in the order 1, 2, 3, 3*, 2*, and 1* as indicated in 
Figure [4]) , the image curve in the Y-plane encircles the origin exactly once in the positive sense for t > 
and exactly three times in the positive sense for t < 0. See Figure [3] This gives the count of the zeros in 
the domain of analyticity. The proof is complete upon using continuity of Y with respect to t and Schwarz 
symmetry. □ 

It is not difficult to obtain the asymptotic expansion of Y{z; g, q*) as z — > oo with 3{z} > 0. The easiest 

r continue 

$'(s)ds 



way to do this is to analytically continue the representation (85) to z outside the loop L, by extracting a 
residue: 



(91) Y(z;q,q*) = 



1 



(92) 



2ni J L S(s)(s — z) 
spai 



(*r 



$'(z) Vat 



outside of L with ^s{z} > 0. 



S(z) a 2 

Now we can let z tend to infinity by expanding the Cauchy kernel in a geometric series for L fixed. Thus: 
1 f &(s)ds 

T -l*) = — O — :p ■ — - f - - < / | 

1 

72 



2?™ Jl S( s )( s ~ z ) 



* Jl S(s) J z 

7T 



S+(s) 



o 



z 



z- 



z 



oo. 



Assuming that q is chosen so that at least the first of the equations (79) holds, we therefore find that 
(93) 1 I - 4 '" 1 ^ 



2ni J L S(s)(s — z) 



(*)* = ^(2(q + q*)t + ax-t)- + 
or z 
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The final step is to recall the expansions (54) and (108) to find 



(94) 



S(z) 



2ni 

z 



0(z 



-3/2\ 



oo, 3{z} > 0. 



Using (|93J) with (|94j) in (|91j) then gives 

4 



(95) Y(z;q,q*) = ^- + 



{2{q + q*)t + ax-t)+2iri 



O 



y3/2 



-> oo, 9{z} > 0. 



4.3. The endpoint 5 as a function of x and t 

satisfied by the endpoint q of the arc B, given (x, t) e 



Now we return to the conditions (79) supposed to be 
2 . First we dispense with the special case of t = 0. 



Proposition 2. The equations (79) hold for t = and x < x c ^/or any arc B connecting q to zp lying in 
the quadrant given by the inequalities 3{z} > and !R{z} < zp ) with q = 3(2;) and q* = 3(2)* being the two 
branches of the complex part of the turning point curve. 

Proof. Set t = and assume that q = 3(x). Whatever the arc B that is the path of integration in both 
integrals in ( |79[ ) actually is, it is homotopic (with orientation preserved) to the image of the map s = i(y) 
as y decreases from x to —00, and the homotopy avoids the line of phantom poles 5i{s} = zp. Since $'(s) is 
analytic for 3{s} > and since the boundary value S+(s) is analytic for s € B, we may deform the path of 
integration from B to the image of s = 3(2/). In the remainder of this proof, we suppose that this deformation 
is done, and that S{z) has been redefined with its branch cut as this new contour and its Schwarz reflection 
in the lower half-plane. 

Therefore, for p — or p — 1, we may parametrize the contour integrals appearing in ( 79 ) by —00 < y < x: 



(96) 



I P (h(x),i{x)*) 




<J>'(s)s p ds 



is 



&(i(v))i(y) p i'(v) dy 



(To get the sign correct it is important to take into account that S+(s) denotes the boundary value taken on 
the integration contour from the left according to the indicated direction of orientation; since S changes sign 
across its branch cuts the limits of integration in the s-integral can therefore be exchanged without changing 
the sign of the integral.) We will prove the proposition by showing that 



(97) 



2nx 



and 



lMx),i(x)*) = ~ M{q}. 

a 



At this point (having used analyticity of $(s) apparent from the formula (48) obtained from Stirling 
asymptotics to deform contours) it becomes more convenient to use the alternate formula (61 1 for ^(s) that 
is obtained from WKB theory. The WKB formula (61 ) is initially defined for real z for which there exist two 



turning points (in the present case this is the interval zl < z < 0) but for the particular initial data under 



consideration it has an analytic continuation to 3{z} > that agrees exactly with the Stirling formula (48) 



By differentiation of 3>(s) as given by (61), we have 

1 duj 1 r x -( s ) 

(98) <Z>'(s) = - — (s)x_( S ) + - 

a as a 



du} 
ds 



(v;s) 



duj- 
ds 



is) 



dv, zl < 8 < 0, 



where u> and lo_ are defined in (56)-(58), and where X-(s) is the real function defined for zl < s < by the 
portion of the turning point curve illustrated in red in Figure [3] Now, x_ (s) admits analytic continuation 
into the complex s-plane from its interval of definition, and since it has a simple critical point at s = z c (a 
local minimum), there is a curve passing vertically through the point s = z c in the complex s-plane along 
which the analytic continuation of X-(s) is real and decreasing away from s = z c . Clearly, the portion of 
this curve in the upper half s-plane coincides with the contour s = i(y) for — 00 < y < x c . Therefore, when 
s = i(y), we have x_(s) = y, where x_(s) denotes the analytic continuation of the similarly-named function 
from the interval zl < s < 0. Using this information and substituting the analytic continuation of (98 1 along 
the path s = j(y) for —00 < y < x < x c into (96) we obtain 



(99) i-pCaO*),^)*) 



ydy 



1 



, S+Uy)) 

where we have also used the fact that dw_(s)/ds 



a 



dui 



dv dy 



4 (independently of the functions p(-) and u(-)). 
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(100) 



We simplify the first term on the right-hand side of ( 99 ), writing it in terms of a double integral as follows: 

}(y) p d'(y) ydy = fx r i{yYi'{y) dy 4 
a 7_™ S+(t(v)) a 



S+(i(y)) 



i(y) p i'{y) 



Exchanging the order of integration in the double integral gives 



(101) 



i(y) p i'(y) 

s+m) 



ydy 



Ax 
a 



i{y) p i'{y) 
s+{i{y)) 



dy 



Reparametrizing the integrals over y by the map s = this becomes 



(102) 



i{y) p i'{y) 
s+{h{y)) 



ydy 



Ax 
a 



q=i{x) S Pds 



dv dy. 



dy dv. 



S + (s) 



dv. 



Now we simplify the second term on the right-hand side of ( 99 ) , by first exchanging the order of integration 
and then reparametrizing the inner integral by s — i{y): 



(103) 



i(y) p i'{y) 



doj 
ds 



dv dy 



— OO J V 



i{y) p i\y) 
S+{i{y)) 



} (v) S+\ S ) 



ds 
duj 
da 



v; s) - 4 



dy dv 
ds dv. 



Combining these results shows that 



(104) I p ( i (x), 3 {x)*)=H 




q=i{x) sP_ds 



i(v) ds S+(s) 



q=i{x) sPds 



S + (s) 



dv 



(The integral over v cannot be broken up without sacrificing convergence.) Here the paths of integration 
in the s-integrals all lie along the image of the mapping s — l(y), —oo < y < x c . Since S(s*) = S(s)* and 
since for v € K we have ui(v,s*) = u>(v,s)* and hence du>/ds(v, s*) = duj/ds(v, s)* where u>(v,s) is defined 
for complex s by analytic continuation from the interval zl < s < 0, we easily obtain 



(105) I p (3(x),i{x) 

where 
(106) 



2x T . , 1 
^ Jp{x)+ 2^ 



— (v;s) 

i(v) aS 



J p {x) 



S + (s) 
i{x) sPds 



l{v) ' dw , ,sPds 
v;s)- 



3(*)* 



ds 



S + (s) 



4J p (x) 



dv. 



s + ( s y 



Here the paths of integration all lie along the image of the mapping s — $(y), —oo < y < x c and its Schwarz 
reflection (which in particular contains the branch cut of the re-defined function S(s)). 

Now we evaluate all of the integrals over s by contour integration arguments (in particular, this will prove 
that the integrand of the w-integral in ( 105 ) vanishes identically). Since S changes sign across its branch cut, 

s p ds 



and since the latter connects $(x)* to $(x), we have 
(107) J p (x) = 



S(s) 



where the path of integration is a closed, positively-oriented loop that encloses the branch cut of S. We may 
now calculate J p (x) in terms of residues at s — oo. Since 



(108) 



1 



2 

s z 



o 



1 



1 _ 1 

5(s) ~ s 
we see easily that 

(109) J Q (x) = -iTT and J^x) = -iir$t{q}. 

Also, by implicit differentiation of the identity ui(v; s) 2 = 16a 2 spo(v) + (As 



1 + auo(v)) 2 we easily obtain 



(110) 



dcu 

as 



(v;s) 



8a 2 p (v) + 4(4s - 1 + au (v)) 



uj(v;s) 



2(1 



By definition of the complex part of the turning point curve, the roots of the quadratic s) 2 are exactly 
s — %(v) and s — 3(«)* for v € K with v < x c ; thus lu(v; s) 2 = 16(s — }{v))(s — 3(v)*). To determine the 
function uj(v; s) as a function of s for fixed v < x c (that is, to determine the proper way to take the square 
root of u>(v;s) 2 ), we proceed as follows. By definition, oj(v;s) is real and has the sign cr_ for Zi, < s < 0, 
because v lies to the left of both real turning points x±(s) in this case. The value of <r_ is determined from 
(56) under the assumption that z lies in the interval in which there exist two real turning points: z^ < z < 0. 
But ( 56 ) can be written in the form 

(111) a_ = sgn(4(z - z P )) 

so since zp > we have cr_ = —1. Now since v < x c , the roots of uj(v;s) 2 form a complex-conjugate pair 
in the complex s-plane, and as uj(v] s) is a negative real function of s for zl < s < while zp > 0, we 
may construct the analytic continuation from the interval s £ (zl,0) to the complement of a branch cut 
connecting the two roots i(v) and %{v)* coinciding with a sub-arc of the branch cut for S(s), and normalized 
so that <jj(v; s) = 4s + 0(1) as s —¥ oo. With the branch cut for u>(v; s) defined in this way, we can write 
u)(v; s)S+(s) = 4S+(v;s) where S(v; s) :— ui(v; s)S(s) / A = s 2 + O(s) as s — > oo, and S(s;v) is analytic 
except for two cuts: one connecting 3(11) to q — $(x) and lying along the contour s — 3(2/) for —00 < y < x c , 
and the other being the Schwarz reflection of the first one. At last we are in a position to evaluate the 
integral involving duj/ds(v; s) by residues; we first use the fact that S(v; s) changes sign across the contours 
of integration to write 

r»( v )* a,., cPJz l f 8a 2 po(v) + 4(4s - I + au (v)) 



(112) 



>{x) duj, ,sPds 
(v;s)- 



ds 



S+(s) 



dtu s p ds 



s p ds 



/,(*)• «"» S + (s) 2 J AS(v;s) 

where the path of integration is a closed, positively-oriented loop that encircles both of the branch cuts of 
S(v; s). Now since 

®U(v)} 



1 



(113) 

S(v;s) s 

a residue calculation shows that, for p 

f» (s > duj 



1 



${q} 



oil 



0. 



(114) 
and, for p = 1, 



(115) 



sds 
OS b+(s) 



ds 
iW* 



(w;s) 



S + (s) 



9cj 



ds iv > s) S + (s) 



00, 



= — 47ri 



9w .sds 
ds {v ' s) S^s) 



= —Am 



n{q} + n{ } (v)} + -a 2 p (v) 



auo(v) — 1 
4 



-47ri3fJ{gr}, 



where on the second line we have used (78). Using (109), (114), and (115) in (105) we arrive at (97), so the 
proof is complete. □ 



Now we show that the solution of the equations (79) can be continued to nearby (x,t) under certain 
conditions involving the function Y(z;q,q*). 



Proposition 3. Suppose that q is a solution of the equations (79) at x = xq and t = t with ^{q} > and 
for which Y(q;q,q*) 7^ 0. Then the equations (79) have a unique smooth solution q(x,t) valid for (x,t) near 
(x ,t ) for which q(x ,t ) = q. 

Proof. We construct the Jacobian determinant of M and M± with respect to q and q* (here viewed as 
independent complex variables). We first rewrite Mq and Mi in the form 



2ir i 
M (q,q*,x,t) = —(2(q + q*)t + ax- t) - -r 



(116) 



$'(s) ds 



M 1 {q, q*,x,t) = —((3q 2 + 2qq* + 3q* 2 )t + (q + q*)(ax - *)) - 



$'(s) ds 

§'(s)sds 



&(s)sds 



where L is a negatively-oriented contour beginning and ending at zp and encircling the contour arc B exactly 
once, and L* denotes the Schwarz reflection of this contour in the lower half-plane (with orientation induced 
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from L by reflection, that is, L* is positively-oriented). Now holding L fixed, the only dependence in the 
integrals on q and q* enters parametrically through the function S(s). Since by direct calculation, 

5 1 1 , d 1 



(117) 
we find that 

(118) 



dqS(s) 2S(s)(s-q) 



and 



1 



S(s) 2S(s)(s-q*) 



dM _ 47rt i 
dq a 2 8 
8M _ 4?rf i 
~dq*~ ~ ~a^ ~ 8 



&{s)ds f &(s)ds 
l S(s)(s- q ) ~ f L , S(s)(s-q) 

$'(s)ds f $'(s)ds 
l S(s)(s-q*) ~ f L , S(s)(s-q*) 



Comparing with the original formula (85 ) for Y(z; q, q*), we see that 
(119) 



dM tt dM 7T , 

— = -Y(q;q,q) and — = -Y(q ; q, q ). 



Similarly, for partial derivatives of Mi we obtain 
(120) 



dMi tt r „ „ . « 

= -3 6i 9 + 2tg* + oau - t - - 
cc 1 8 



^ = ^ [6tg* + 2tq + ax - t] - 
dq* a z 8 

The partial derivatives of Mi can be rewritten as 

= —= [otq + 2tq +ax — t\ — — 



ft' 



$'(s)sds 
l 5(s)(s - 9) 

$'(s)sds 
L S{s)(s - q*) 

&(s) ds 
L S( S )( S -q) ' 



$'(s)sds 
u* S{s)(s - q) 
<&'(s)sds 

L* S( S )(s-q*) 

$'(s) ds 

L> S(s)(s-q) 



(121) 



&(s)ds f <f>'(s)ds 
l S(s) ~ f L , S{s) 
iq 



dMi it .„ * „ 

[6tq* + 2tq + ax - t] 
or 8 



dq* 



&(s)ds f &(s)ds 
l S(s)( S ~q*) ~ f L . S(s)(s-q*) 



$'(s) ds 



$'(s) ds 



S(s) 



Setting x = xq and t = to, we use the fact that q is a solution of Mq = to eliminate the terms on the 
second line in each case: 



dMi 



(122) 



dq 



dMi 



47rto iq 



Oo,*o) 



$'(s) ds 



dM 



dq* 



dq 



(a:o,to) 

iq* 



<f>'(s)ds 



$'(s)ds / $'(s)ds 



dM 



Oo,*o) 



Therefore, the Jacobian determinant is 

'dMi dM dM dMi 



J r (q,g*,x ,t ) 



(123) 



dq dq* 



dq dq* 



(2:0 ,*o) 
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8 



Y(q;q,q*)Y(q*;q,q*)(q-q*) 



-\Y(q;q,q*)\ 2 Q{q}. 



Since this is nonzero by assumption, the solution to the simultaneous equations M = Mi = can be 
continued from q to q{x,i) by the Implicit Function Theorem. □ 



Corollary 1. Let x < x c . Then, for \t\ sufficiently small, there exists a unique solution q = q(x,t) of the 
equations (79) satisfying q(x,0) — $(x), where 3(2:) is the complex part of the turning point curve. For t > 
the solution can be continued (even for x > x c ) until ^s{q(x,t)} becomes zero. 

Proof. From Proposition [2] we have that the equations ( |79"| ) are satisfied by q = 3(2;) when t = and x < x c , 
and we note that 3(3(2;)} > for x < x c . By Proposition [I] the only possible zeros of Y(z;q, q*) in the 
open upper half-plane are far from q when \t\ is small, which implies that Y ($(x); $(x) , $(x)*) 7^ 0. Moreover, 
if t > then Y(q;q, q*) can only become zero if = 0. Therefore, by Proposition [3] the Corollary is 

proved. □ 

It follows that for t > 0, the only obstruction to smooth continuation of the endpoint function q(x, t) from 
its values q(x, 0) = 3(2;) for x < x c may be a collision of q with the real axis. We now completely resolve the 
scope of the continuation. 

Proposition 4. There exists a well-defined smooth curve x = x c (t) for t > satisfying x c (0) = x c , such 
that the solution of (79) with q(x,0) = 3(2:) for x < x c can be uniquely continued to the domain t > and 
x < x c (t). Also, ^s{q(x, t)} > whenever x < x c (t), but ^s{q(x, t)} 4- as x f x c (t). 

Proof. First, we argue that it is impossible for ^s{q} to vanish for bounded (2;, t) unless q tends to a real 
point z in the interval (zl,0). Consider the combination of Mq and M\ given by: 

F(q, q*;x,t) :=M 1 (q, q*; x, t) - \{q + q*)M (q, q* ; x, t) 

(124) irt 1 

a z 2 

Clearly, we must have F(q,q*;x,t) = for any solution of the two equations Mq = Mi = 0. Now for t 
bounded the first term on the second line above obviously converges to zero as | 0. By definition of 

Ip{q, q*) we have 



(125) h(q, q*) - -(q + q*)I {q, q*) = Z j ^ &(s) £ ' ds 

A careful dominated convergence argument (in which the integration from q is replaced locally by a loop 
around s = q that is deformed to the real axis in a neighborhood of z £ R, a contour on which (s — \{q + 
q*))/S(s) converges pointwise to 1 and is uniformly bounded) shows that 



h(q,q*) - ^(q + q*)io{q,q* 



£ P $' + ( s )rf s | = -3{$+(z)} 



(126) lim 

9-S.26R 

O{<?}>0 

because 9{$ + (zp)} = according to (53). But again according to (53) we see that 9{<i> + (z)} = only for 
z — zp and for z € [zl>0]- Therefore, for bounded t we can only have 4- if q — » z e \z\, , 0] or if q — > zp . 
But a similar argument also shows that Mo(q, q*; x, t) — ¥ 00 if q — > zl, q — > 0, or if q — >• zp with > 

and (x,t) bounded, so q can only tend to a point z G (zl,0) for bounded (x, t). 

Now fix z £ (z\j, 0). Then F(q, q*; x, t) — > as q — > z with > 0, so the two conditions Mq = Mi = 

degenerate to a single condition: M§(z;x,t) — 0, where 

2tt 

(127) Mo(z;x,t) := (4zt + ax ~ t) + Iq(z), z h < z < 0, 



where 



(128) 



where on the second line the path of integration from zl to zp lies in the upper half s-plane. Since $'(z) — > 
+00 as 2 4- z h while $'(z) — > -co as z f 0, the expression on the first line shows that Iq{z) — > —00 as z 4- 2l, 
while that on the second line shows that — ^ as z t 0- Both expressions are real- valued, finite, and 
equivalent for z G (zl, 0). 
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For each z £ (zl,0), the condition Mq(z;s,<) = defines a straight line in the (x, i)-plane. We will now 
show that the union of these lines is an unbounded region of the form x > x c (t) for all t > 0, where x c (i) 
is a smooth function of i with x c (0) — x c . More precisely, through each point (x,t) with x > x c (t) and 
t > there will pass exactly two of these lines corresponding to distinct values of z £ (z L , 0), and the curve 
x = x c (t) will appear as the caustic formed by the intersection of infinitesimally neighboring lines. Given 
t > fixed, the possible x- values lying on the lines are given by 

1 — 4z a 

(129) x= t--]J(4 z L < z < 0. 

Of Z7T 

The first term is obviously bounded for t fixed and z £ (zl,0), but the second term tends to +oo as z 1 zl 
or z j" 0. This shows that for t > fixed the region occupied by the lines is exactly 



(130) 



x > x c (t) 



mm 

z L <z<0 



1 -4z 



-t 



2tt 



< +oo. 



To determine x c (t) we show that for each t > 0, the function to be minimized has exactly one critical point 
in (zl,0), which is then (by asymptotics as z 1 zl and as z f 0) necessarily the unique minimizer. Indeed, 



(131) 



d_ 
dz 



1-42 



i - 



2tt 



^(z;x,t) 



2?r <9z 



~Y c (z;z), 



where Y c (w; z) is by definition for zl < w < the limiting value of the function Y(w;; g, g*) as g — > z £ (zl, 0) 
with ^{q} > 0: 

$'(w)-$'(z) 2 ( 



(132) r c (w;z) 



16t 

A 2 



$' + (s) 



zl < w < 0, zl < z < 0. 



(s — z)(s — w) 

It is obvious that Y c (w; z) is symmetric in its arguments: Y c (w; z) — Y c (z; w), and for each z £ (zl, 0) there 
is, by Proposition [T] a unique simple zero w — w(z) £ (zl,0) of Y c (w;z). The graph w — w(z) must be 
symmetric with respect to reflection through the diagonal w — z because otherwise there would necessarily 
be multiple roots w £ (zl,0) of Y c (w;z) for at least some z in the interval (zl,0). Therefore, this graph 
either coincides with the diagonal line w = z or it connects (zl,0) with (0, zl) and crosses the diagonal 
exactly once orthogonally. But the case w(z) = z can be ruled out easily by symmetry of Y c (w; z) and the 
fact that Y c (w; z) f +oo as w i zl while Y c (w; z) 1 — oo as w t 0. Therefore, Y c (z; z) has exactly one simple 
root in the interval (zl, 0) for all t > —T, and this establishes the uniqueness of the critical point and hence 
of the minimizer. Let us denote this critical point by z = z c (t), and the corresponding value of x defined for 
z = z c (t) by (129) as x — x c (t). By consideration of a different type of configuration for the (/-function in 
which the branch points of the square root S(z) are both real, and examination of the limit in which these 
two points coalesce (see [9]), it can be shown that in fact z c (0) = z c £ (zl, 0), where z c is the specific value 
at which the complex branches $(x) and %(x)* of the turning point curve become become real (see Figure p|. 
It then follows from Proposition [2] that x c (0) = x c . □ 

The curve x = x c (t) is plotted for the case of a = S = 1 and /i = 2 in Figure [5j In Figure [6] we plot the 
differences between x c (t) for various parameter values and the specific function x c (t) shown in Figure [5] In 
particular, these plots clearly show that x c (t) is not (in general at least) a linear function of t. 

Finally, we show that the complex- valued function q(x, t) defined for alH > and x < x c (t) satisfies certain 
canonical partial differential equations (Whitham equations) in the interior of the domain of definition. 



Proposition 5. Suppose the equations (79) admit a solution q — q(x,t) differentiable with respect to x and 
t . Then 



(133) 



dq 
dt 



-[1-3,-^=0 
a ox 



dq* 
dt 



Proof. This follows from the relations Mo(q(x,t),q*(x,t),x,t) = Mi(q(x,t),q*(x,t),x,t) = by implicit 
differentiation. Indeed, since 

9M 2tt 8Mq 2tt 

(134) -dx- = ^ and ^r = ^ [2((?+(z) - 1] 
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Figure 5. The curve x c (t) for a = S = 1 and [i = 2. The four marked points in the domain 
x < x c (t) correspond to the (x, i)-values in the four indicated panels of Figure [7] below. 
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Figure 6. Left: for a = S = 1, curves x = x c (t) for various values of \i compared with the 
curve in Figure [5] Right: for a = 1 and /i = 2, curves x = x c (t) for various values of S 
compared with the curve in Figure [5] 



while 
(135) 



9Ml % + and 9 ^ = ^[3<? + 2 qq * + 3q* 



dx a dt 
it follows immediately from (119) and (122) given in the proof of Proposition [3] that 

(136) & = 

and also that 

dg _ 4 . , , <V 4 

9t a 2 F(g;g,<7*) 



(137) 



a?F(<7;g,<7*) 9x 



ar(g*;g,g* 



4 <9g* 

1 - 3g - Q* and — = [1 - 3<f - g . 

9i a 2 y(g*;g,g*) 



It is then obvious that ( 133 1 holds true if Y(q; q, q*) and Y(q*; q, q*) are nonzero, a condition for the existence 
of a smooth solution q(x, t). □ 



2.i 



The equations ( 133 ) are actually equivalent to the dispersionless MNLS system (|9| in the locally subsonic 



case that Q < 0, as can be seen by means of a straightforward calculation using the substitutions 

(138) p =±(%{ q V*}f and U= I (1 _ 4 | 9 |). 

or a 

This shows that the variables q and q* are complex Riemann invariants for the elliptic dispersionless MNLS 

system. 

4.4. Choice of the arc B. Supposing that t > and x < x c (t), we determine the endpoint q = q{x,t) 



with 5{g(x,i)} > as explained in £4.3 and hence obtain as explained in [4.2 the functions g(z) and h(z) 
(parametrized by (x, t) of course). These functions are actually only well-defined once we specify a particular 
contour arc B in the upper half-plane connecting z = q with z = zp . We now describe how B is to be chosen. 



Proposition 6. Let t > and x < x c (t), and let q = q(x,t) be determined as explained in {4-3 Then there 
is a unique choice of the arc B connecting z — q with z — zp for which 9{/i(z)} can be extended by continuity 
to the domain z G C \ (-co, Zi] U [0, +oo). 

Proof. The key observation is that as the imaginary part of an analytic function, 3{/i(z)} is a harmonic 
function in the domain C \ ((-co, zi] U [0, +oo) U B U B*). Generally there is a jump discontinuity across 
BUS* that we wish to avoid by choice of B. 

Firstly, regardless of how B is chosen, we can easily see that 3{/i(z)} = for z G [zl,0], for z = zp, 
for z = q, and for z = q* . Indeed, by Schwarz symmetry of h and the fact that (zl,0) is in the domain 
of analyticity of h, it follows immediately that Q{h(z)} = for zl < z < 0. But since S(z) is Schwarz 
symmetric and analytic on the real axis except at its jump discontinuity point z — Zp, and since according 
to (|90|, the boundary values 3{Y± (z; q, q*)} taken by 3{F(z; q, q*)} on the real axis are bounded, it follows 



from ( |89[ ) that the corresponding boundary values of 3{/i(z)} are real differentiable functions of real z 
except at the point z = zp (which is nonetheless a point of Lipschitz continuity of the boundary values 
of 3{/i(z)} for z G M). In particular this implies that the limit points z = zl and z = of the open 
interval (zl,0) are also points where 3{/i(z)} = 0. The fact that the well-defined value 3{/i(zp)} vanishes 



is a consequence of the Fundamental Theorem of Calculus, the formula (89), and the formula (90); since 
S(z) — —y/{z — 3?{g}) 2 + 5{g} 2 holds for z < zp, we integrate along the top edge of the branch cut for 
Y(z; q, q*) to obtain 







(139) = \ \/{s ~ ®ti}) 2 + %{q} 2 %{Y + (s; q, q*)} ds 



= -2tt / - + - / ds. 

Jo vV - 16s a Jz+ 

By direct evaluation of these integrals and the use of the definitions of z + and zp we obtain the claimed 
result that 9{/i(zp)} = 0. By Schwarz symmetry of Y(z;q, q*) the same result holds had we integrated 
instead along the lower edge of the branch cut for Y(z; q, q*). With this information we can express 3{/i(g)} 
as an integral again using the Fundamental Theorem of Calculus. Indeed, integrating along the left edge of 
the branch cut B in the direction of its orientation from z = q to z = zp, we get (using 3{ft,(zp)} = 0) 

(140) %{h(q)} = ^ S + (s)Y(s; q, q*) ds 



Now using the formula (88) and the fact that g(z) is Schwarz-symmetric and analytic for z G C \ (B U B*) 
we have 

(141) = \^ {J b (9' + ( s ) - 9'-(s)) ds} = \jg'{s) ds=~ jp{s) ds, 

where the closed contour of integration is a large positively-oriented circle of arbitrarily large radius. It 



therefore follows from the fact that q satisfies the conditions (79) that p(s) = 0(s~ 2 ) as s oo and hence 
= as claimed. The fact that 3{/i(<7*)} = then follows by Schwarz symmetry of h. 
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Still regarding the arc B connecting z = q with z = zp as arbitrary, we note that the zero level set of 
3{h(z)} defined by Z := {z e C\ ((-oo, z L ] U [0, +00) UBU B*), ^{h{z)} = 0} actually extends by continuity 
to the excluded branch cut BU B* and moreover is independent of B. This follows from the fact that given 
Sy{/i(q)} = we may express ^s{h(z)} in terms of an integral from s = q to s = z, and because s = q is the 
square-root branch point of S(s), a change of choice of branch for S simply amounts to a change of sign of 
the integral and hence of ^{h(z)}. This clearly leaves Z invariant. 

We will now characterize Z completely, and prove in particular that Z contains a smooth arc connecting 
z = q with z — zp. When the proof is finished we will select this arc to be the arc B and show that this 
choice makes ( =s{h(z)} continuous. We firstly characterize the real points of Z. We already know that Z 
contains the real interval [zl,0] and the real point zp. Moreover, it is easy to see that these points exhaust 
ZflR. Indeed, since S(z) is bounded away from zero for real z^zp and sgn(5(z)) = sgn(z — zp), it follows 
from (89) and (90) (see also Figure [4| that is strictly increasing for z + < z < zp and strictly 



decreasing for z < zl, for < z < z + , and for z > zp. 

Next we consider the points of Z in the upper half-plane near the real axis. According to Proposition [l] 
for t > there exists a unique simple zero z = £ of Y(z; q, q*) in the interval (zl, 0), and since S(z) is nonzero 
near this point, z — £ is obviously a simple saddle point of 3{/i(z)} and therefore there is a unique branch of 
Z emanating from this point with a vertical tangent into the upper half-plane. Also, since the real derivative 
of the boundary value 3{/i + (z)} has a jump discontinuity at z = zp with opposite nonzero left and right 
limits, there is a unique branch of Z emanating transversely to the real axis into the upper half-plane from 
z = zp (here the tangent is not necessarily vertical, however). The points z = £ <E (zl,0) and z = zp > 
are the only real limit points of the part of Z in the open upper half-plane. 

Next, since z = q is a simple root of S(z) 2 , a local analysis of S(z)Y(z; q, q*) near z = q using the fact (see 
Proposition [I]) that Y(q; q,q* ) 7^ shows that there exist exactly three arcs of Z emanating from z = q € Z 
separated by angles of 2tt/3. 

Let us now analyze Z assuming that \z\ is large with 9{z} > 0. Since S(s) = s[l— \{q + q*)s~ 1 +0(s~ 2 )] 



as s —> 00, we recall (89) and (95) to obtain 



9{M*)} = -^{ [* S(s)Y(s;q,q*)ds 
(142) " U " 



'1 

At , ( 2 



-z 2 - — (ax - i) + vri I z + <3(z 1/2 H , z -> 00, 3{z} > 0. 



For t = 0, the part of the level set Z in the distant upper half-plane is clearly a smooth curve asymptotic 
to the straight line through the points z = and z = ~2x/a + iri as z — > 00. For t however, Z is 
asymptotic to leading order to the union of the real and imaginary axes, and the computation of higher order 
corrections is necessary to determine whether the horizontal asymptotes actually correspond to branches of 
Z in the upper half-plane. Representing a branch of Z in the form z = \z\e z ^, 4> = ( t ) (\ z \)i we substitute into 
( |142| ) and set Q{h(z)} = to obtain the relation 

(143) - 4sm(20)- -^—(ax-t)sm^)- ^ cos(<?()) + 0(|z|- 3 / 2 ) = 0, \z\ -> 00. 

a z a z \z\ \z\ 

Setting |z| = 00 gives sin(2</>) = 0, all roots = <fio of which are simple. We may therefore apply the Implicit 
Function Theorem to continue these roots to finite |z|, yielding the asymptotic expansion 

(144) m) = ^0 - { Zn ] siI S\ - ^t^tSt + °(N- 3/2 >' 

4t\z\ cos(20 o ) ot\z\ cos(20o) 

where 4>q is a root of sin(20o) = 0. Since we are only concerned with the upper half z-plane, we need to 
consider only the angles (f>o — 0, ^7r, n. For 0o — the asymptotic formula predicts a slight deformation of 
the angle in a direction depending on x and t, but regardless this branch remains in the upper half-plane due 
to the dominant constant term in the angle <f>. For O = and 4>q — n, we have sin(0 o ) = 0, and the surviving 



term at order |z| 1 in the asymptotic expansion (144) indicates that for t > we have \ ^ir — 4>(\z\)\ > for 



large \z\, showing that these solutions do not correspond to branches of Z. Therefore, we conclude that for 
t > the part of Z in the distant upper half plane consists of exactly one smooth curve tending to infinity 
asymptotic to a non-horizontal straight line. (If t < then the same analysis shows that when (f>o — or 
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0o = 7r we have instead \\n — 4>(\z\)\ < for large \z\ and there are then three distinct branches of Z in the 
distant upper half-plane.) 

So, we know that in the open upper half-plane there exist branches of Z emanating from the real axis only 
from z = £ G (zl, 0) and from z — zp. There also exist exactly three branches emanating from z = q, and 
(for t > 0) there is exactly one unbounded branch in the upper half-plane. Since Q{h(z)} is harmonic for 
> except along the branch cut B across which 9{/i(z)} simply changes sign, the Maximum Principle 
prohibits the existence of any isolated components of Z in the open upper half-plane, with the possible 
exception of curves enclosing z = q. But the latter are also easily excluded from Z by local arguments. 
Therefore, the part of Z in the upper half-plane consists exactly of the six branches described above, and 
these must be matched with each other. The Maximum Principle also prevents any of the arcs emanating 
from z = q from coinciding with one another. Finally, since Proposition [T] rules out any zeros of Y(z;q,q*) 
in the open upper half-plane for t > 0, there are no saddle points of $s{h(z)} for > and this means 

that the six arcs must be matched pairwise, and no crossings are allowed. Since none of the three branches 
emanating from z — q can be matched with each other, we conclude that exactly one of these three branches 
is connected to each of the three distinct points z — £, z — Zp, and z = oo. As the zero level set is 
Schwarz-symmetric, this concludes the complete description of Z in the case t > 0. 

Now suppose that the branch cut B is taken to coincide with the branch (a smooth arc, actually) of 
Z connecting z — q with z — zp. It then follows that 3{/i(z)} can be continuously extended to z € 



B U B* . Indeed, since 3{ft,(<7)} = 0, from (89) and the Fundamental Theorem of Calculus, the mismatch 



9{/i + (z)} — 5{/i_(z)} of boundary values taken by ^s{h(z)} on B is 

3{Mz)}-Q{M*)} = -is<; / s+is)Y(»:<i.<r)<i»- I .s'_ ( .v.r ( .s: </.,/. </.s 

(145) 



l q J q 

. sGB seB 



^{j S+(s)Y(s;q,q*)ds> . : G B. 
.seB 

But this is the same as 23{/i + (z)}, which vanishes because z £ B C Z. □ 

From now on we assume that for all t > and x < x c (t), the arc B is chosen exactly as explained above, 
making 3{/i(z)} continuous across B U B* . The zero level set Z clearly divides the open upper half-plane 
into three disjoint open regions. We denote by V (the valley of Q{h(z)} in the upper half-plane) the union of 
the two open regions that are separated by i3, and we denote by M. (the mountain of Q{h(z)} in the upper 



half-plane) the remaining open region. According to ( 142 ) and the fact that Z contains all of the points for 



which 3{ft,(z)} = 0, we have the following strict inequalities for > 0: 

(146) ^{h(z)} > 0, z e M 
and 

(147) 3{/i(z)} < 0, zeV. 

The regions A4 and V are illustrated for several choices (x,t) in Figure [7] 



4.5. Key properties of g(z) and h{z). With g(z) and h(z) now completely determined as analytic functions 
depending smoothly on parameters (x,t) with t > and x < x c (t), we now explain all of the important 
consequences of our finished construction. 

Proposition 7. The function g is analytic for z £ C \ (BUB*), satisfies the Schwarz symmetry condition 
g(z*) — g(z)* , and also g(0) = while g(oo) is well-defined and finite. The function h is analytic for 
z G C \ ((— oo, zl] U [0, +oo) U BU B*), satisfies the Schwarz symmetry condition h(z*) = h(z)* , and 3{/i(z)} 
is continuous for Q{z} =/= 0. There is a real constant k — n(x,t) such that 



(148) 



h+(z) + h-{z) = k, zeBuB*. 
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Figure 7. Four numerically generated sign charts for Q{h(z)} indicating the regions M. 
and V, all for parameter values a = 1, S = 1, and /i = 2. Upper left (UL): (x,t) — (—1, 1). 
Upper right (UR): (x,t) = (0.25,0.5). Lower left (LL): (x,t) = (-1,0.5). Lower right (LR): 
(x,t) = (—0.5,0). The red curves are branch cuts for h(z), the blue curves are zero level 
curves of 3{/i(z)}, and the purple curve is B U B* , which is both a branch cut of h(z) 
(although ^s{h(z)} is continuous there) and a zero level curve of ^s{h(z)}. 



Also, the strict inequalities (146) -(147) hold for > 7 with the only excluded points in the upper half- 

plane satisfying 3{/i(z)} = 0. From (142) we have the asymptotic condition 



(149) 



%{h(z)} = - 



2tt 

/"V —* 



St 



%{z} 



^{z} + 0(\^{z}\^ 2 ), 5R{z} -> -ex), < = O(l). 



Finally, h{z) is continuous at z — q, and 

(150) h\zf 



(z - q)v(z) 



where v is analytic and non-vanishing in a neighborhood of z = q. 



Proof. It only remains to prove (148]) and (150). But, these both follow from the representation (89). 
Indeed, taking into account that S{z) changes sign across B while Y(z;q,q*) is analytic in a neighborhood 
of B, we easily see that the sum of boundary values h'Az) + h'_(z) vanishes identically for z G B. Therefore, 
h + (z) + h-(z) is constant along B, and since B is part of the zero level set Z of ^s{h(z)}, it is obvious that 
the constant value taken by h + (z) + h-(z) for z g B is real. It then follows by Schwarz symmetry of h(z) 
that the same identity holds for ze8*. This proves (148). The formula (150) also follows directly from (89) 
using Proposition [l] □ 

With the parametric dependence of h on (x, t) completely determined, we can consider the partial deriva- 



tives of h with respect to x and t. Recalling the definition (84) of h in terms of g, 9, and <&, and the fact that 



^(z) is a function independent of (x,t), we immediately deduce that these partial derivatives are analytic in 
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a larger domain than is h itself, namely for z £ C \ (BU B*). We claim that the partial derivatives are given 
by the following simple explicit formulae: 



(151) 

and 

(152) 



dh 

dx 



1 8k 
2dx 



-S(z), z eC\(BuB*), 



dh 



(2) 



1 dn 
Ydt 



1 



* + «{?}-- zeC\(8u8*). 



Indeed, these formulae exhibit the correct domain of analyticity and capture the correct boundary conditions 

8k 



on BUB* following from ( |148 ) 
dh 



(153) 



and 



zeBUB* 



dx w <9a; dt K ; ' dt y ' dt 1 

Also, since the partial derivatives of g are bounded for z = oo, the principal parts of the partial derivatives 
of h at z = co must agree with those of 0, and by expansion of S(z) for large z it is easy to confirm that the 
above formulae indeed satisfy 



(154) | W 



8e_ 

dx 



(z)+0(l) 



-2+0(1) and 



dh 

~dt 



at 



(z)+0(l) 



t 2 

^z 2 H — -2+0(1), z-^oo. 



Since g(0) = for all (x,t) it follows from (84 1 that the partial derivatives of h must agree with those of 
exactly at z = 0. Imposing these conditions on (151) and ( 152[ ) proves that the partial derivatives of k can 
be expressed explicitly in terms of q(x,t). Indeed, since 5(0) = — \q\, 



dh .. d9 1 dn 



(155) 
and 



dh , , 89 . . 1 8k 
^ °=8t {0) ~8t {0) =2 8t 



a 



1 

2a 



»{«}- 



1 

4a 2 



— = -(l-4 9 ) 
dx a 



| = A (1 - 2K{g})M 



1 

2a^ 



Similar calculations allow us to establish a direct expression of k = K,(x,t) itself. Indeed, since g(z) is 
analytic for z e C \ (S U B*), is well-defined for z = oo, and satisfies g+(z) + g~(z) = 29 (z) + $(2) — k for 
z e BU B*, it follows that it must be given by 

(157) 



S(z) 
2m 



BUB* 

1 

2' 



26>(s) + $(s) - k 
S+(s)(s - z) 



25(2), 

= 6> (z) --K+ — -7T 1 {ax - t) + 



ds 

AtS{z) 



(z + R{q(x,t)}) 



S(z) 
2-ki 



$(s) 



z eC\(Bu8*). 



Enforcing the condition that g(0) = and using the fact that S(0) 
k(x, t): 



bjb* S+{s)(s-z) 

-\q(x,t)\ then gives a formula for 



(158) K(x,t) 



2a 



'" {2ax-t)-^\q{x,t)\(ax-t)-^\q(x,t)\^{q{x > t)}-^^ 



$(s) ds 

BUB* S S+( S ) 



In particular, for t — we have 



(159) 



k(x,0) = -(l-4|g(x,0)|) 



\q(x,0)\ 



$(s) ds 

BUB* s S+i s ) 



where in the last integral we use q = q{x, 0) and the corresponding arc B. This formula can be simplified 
even further by taking (carefully) the limit x f x c in which q(x, 0) — > z c < 0: 



(160) 



lim k(x,0) = — (1 + 4z c 



2$(0)-$(z c ) 



2z r 



3{$ + (s)}ds 
s(s - z c ) 



Note that $(0) is a well-defined real number, and the fact that the integral is convergent at s = follows 
from (53) which shows that 3{$ + (z)} vanishes linearly as z I 0. 
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5. Steepest Descent Analysis of the Riemann-Hilbert Problem 



5.1. Opening of Lenses and Introduction of g(z). We begin by making an explicit modification of the 
matrix N(z) by "opening lenses" according to the diagram in Figure p] The most important features that 




Figure 8. The domains A, £, and 1Z, along with their mirror images in the lower half- 
plane. Note that dC n dlZ coincides with the contour arc B. The solid curves comprise a 
contour denoted S. 

we require of the contour in this figure are: 

• The contour arc C lies entirely within the mountainous region A4 for Q{h(z)}, and is asymptotically 
horizontal as 3?{z} — » — oo. 

• The contour arc B+ lies entirely within the unbounded component of the valley region V for Q{h(z)}. 

• The contour arc B- lies entirely within the bounded component of the valley region V for 

• The parameter 77 is sufficiently small that zp — r\ > z + . 

The modification is defined as follows: 

(161) O(z) := N( z )e t9 ^ x ^/ e 

(162) O(z) := N(z)e ie(j:;:E ' t),T3/e 

(163) O(z) := N(z)e lf,(z; ^ t)CT3/£ 

(164) O(z) := aiO(z*)*a u zeA'UTUK', 

and elsewhere we set O(z) := N(z). Recalling that E{z) is an analytic nonvanishing function for z in the open 
upper half-plane, we see that this defines the matrix O(z) as an analytic function of z in the complement of 
the system of contours illustrated in Figure |] Then, we introduce the function g(z) defined in [j4]by setting 

(165) P(z) := 0(z)e t9{z)a3/e . 

Since g(z) is analytic for z € C \ B U B* , and since B = dC n d7Z, P(z) has the same domain of analyticity 
as does O(z). 
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-£(z)e 4 *( z )/ £ 
1 



^ — i9(z;x,t)<T 3 /e 



z e A, 



-E{z) 



-l e -i*(z)/e 



E{Z) 



C3/2 — i6(z;x,t)cr 3 /e 



ze C, 



E{z)-\ 



»0(*)/e 



-E{z)e l ^ z >l e 
1 



E{zY 3/2 e~ l0{z ' x ' t),j3,t 1 zeTZ, 



The explicit transformation of N(z) into P(z) differs from a more standard choice that is based upon the 
factorization of the jump matrix for N(z) (say, for z < 0) given by 



(166) 



s(zY 



1 

-2i9(z;x,t)/e 



l-\s(z)\ 2 



1 - s ( z )e 2ie( - z ' x ^ /e 
1 



This more standard choice was used, for example, in the analysis of Tovbis, Venakides, and Zhou of the 
semiclassical limit for the focusing nonlinear Schrodinger equation |21j with special initial data for which 
there were no discrete eigenvalues (poles of N(z)) and for which the reflection coefficient analogous to s(z) 
was given in terms of Euler gamma functions. It has the advantage that after opening lenses there remains 
no jump across the real axis for z < zp. However, opening lenses based on the factorization (166) actually 



sets into motion a chain of unfortunate events leading to more difficult analysis and even some ambiguity 
about the scope of the resulting asymptotic formulae. The main problem with using the factorization ( 166 ) 
to open lenses is the presence in the upper half-plane near z — zp of the phantom poles of s(z). These poles 
both constrain the possible location of the contours near z = zp (to avoid introducing new singularities 
into the Riemann-Hilbcrt problem for P(z) jump matrices involving s(z) must not be deformed through 
the phantom poles) and also lead to nonuniformity in the Stirling asymptotics of s(z) near z = zp that 
ruins the approximation E(z) w 1. The authors of |21) dealt with the latter issue by installing an unusual 
local parametrix near the point analogous to z = zp in their problem; the parametrix is difficult to analyze 
because it is not explicit (its existence relies on abstract Fredholm theory) and it depends on e in a way that 
must be carefully understood. By contrast, the reader will see that the approach we use here (which is based 
on extending not s(z) but rather the product (1 + e 2 ^ z ^ t )s{z) — related to the function E{z) by ( |5l] ) - 
into the upper half-plane) removes all difficulties with the phantom poles near z = zp and moreover avoids 
completely the need for any sort of local parametrix near this point. In fact, E(z) is analytic throughout 
the upper half-plane, and it will turn out that in our approach the contribution to the error from the 
neighborhood of z — zp is exponentially small. 

We recall that g(0) — and also that E + (0) = (the latter following from the fact that s(0) = while 
<&(0) is well-defined and finite) and hence it follows that P(z) satisfies the normalization condition 



(167) 



P(0) = 



(the value is the same whether the limit is taken from the upper or lower half-plane). Since 3(00) is well- 
defined, P(z) tends to a diagonal limit as z — > 00 for each fixed e. To describe the jump discontinuities of 
P(z) across the arcs of the jump contour pictured in Figure [8j we suppose that all contour arcs are oriented 
left-to-right (that is, from —00 toward q, q* , or zp — rj, from q or q* toward zp — n, zp, or zp + 77, from zp — n 
toward zp, from zp toward zp + 77, and finally from zp + rj toward +00). 

Then the jump conditions satisfied by P(z) along the real z-axis are as follows: 



(168) 



P + (z) = P_(z) 



1 e 2f{z)/e s f z \ e 2i{e{z;x,t)-g{z))/e 
,2f(z)/e s r z y e -2i(e(z;x,t)-g(z))/t 1 _ e 4/( Z )/e| s ^\|2 



Z < 0, 



(169) P+(z) = r 3 P_(z)?r <T3 



1 e 2f(z)/e s ^ e 2i(9(z-x,t)-g(z))/e 
e 2f(z)/e s f z y e -2i(6{z;x,t)-g(z))/e 1 + e 4/(z)/e| s ( 2 ) |2 

< Z < Zp — 77, 



(170) P+(z) =i CT3 P_(z)r CT3 
where 



A n (z) Al2 r z yi(9(z;x,t)-g(z))/t 

A 21 (z)e- 2t< - e ^ x - t ^9(-))/e A22 ( z ) 



zp — 77 < Z < Zp, 



(171) a,i(z) ■= 1 + |s(z)|2 e s{*+(*)}A 



(172) 



*12(Z) A 21 (Z) . (1 + e2/(z )/ e) | s(z) | e 
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and where Au(z)A 22 {z) - A 12 (z)A 2 i(z) = 1, 

Bu(z) 



(173) P+(z) = i' T3 P_(z)r ,T3 
where 

(174) B 22 {z) 



£ 12 (,z) e 2 »(0(*;*.f)-ff(*))A 

5 21 ( z ) e - 2 «( e (^*)-9(*))A 5 22 (z) 



Zp < Z < Zp + T], 



1 + 



(l + e 2/WA)| s (^)| 



K*+(z)}A 



(175) 



12 (z) = B 21 (zY •= M')^" ~ 1 e*{« + M}A 



and where Bn(z)B 22 (z) — B 12 (z)B 2 i(z) = 1, and finally 

r i _^f z \ e 2i(e{z-x,t)-g(z))/(.- 

(176) P + (z) = f 3 P_( Z )i ff3 [_ s(z) , e _ 2l(fl(z;I , t) _ 9(z))/£ l + | s (z)| 2 J' z>z p+ ? ?- 

These jump conditions may appear complicated, and this is the price to be paid for eschewing the standard 
two-factor factorization ( 166 1 in favor of a more complicated one. However, we will now prove that all of the 
jump conditions for P(z) along M amount to exponentially small jump discontinuities, given the condition 
j] < zp — z + . Indeed, using (42) and (43) and recalling that 9(z;x,t) and g(z) are real for z £ R, we easily 
see that ( 168 1 takes the form 

(177) 



P+(z) =P_(z) (l + 0(e 2/ W/ £ )) , z<0, 



that ( 169 1 takes the form 
(178) 



P+(z) = i CT3 P_(*)r CT3 (l + Ofe^ 6 )) , 



and that (176) takes the form 

(179) P+(z) = i CT3 P_(z)r CT3 fl + 0(e" /(z)/e 



< z < zp — 77, 



z > zp + <q. 



In each case, the jump matrix is a uniformly exponentially small (in e) perturbation of the identity that also 
decays exponentially as z — > ±00. Also, due to the estimate (44), we have that 



(180) 

and 

(181) 



P+(z) = P_(z) (l + OG^/W^)) , -77 < z < 0, 



p+(z) =r 3 P_(^)r CT3 (i+og^/W^)) , o<z<ri 



both hold uniformly for some 77 > and M > 0. Next, since according to (53) we have ^{^^-(z)} = — |/(z) 
for z > z + , we sec that for zp — 77 < z < zp we have 



(182) 



V 7 (1 + e 2 /( z )A)|s(z)| 

(|s(z)|e^ z )/ e - l)(|s(z)|e^ z >/ e - e 2 ^ z )/ e ) 
+ |s(z)|e/WA(l + e 2 /W/ e ) 

= 1 + exponentially small in e, uniformly for zp — 77 < z < zp, 



where we have used (47) and the fact that e 2 ^ z -*/ £ < 1 for z < zp. Similarly, from (47), 

|1- |s(z)| 2 e 2 ^ z )/ e | 



\A 12 (z)\ = \A 21 (z)\ = 



(183) 



< 



2|s(z)|e/W/ e cosh(/(z)/e) 
|1- \s(z)\ 2 e 2 f^/ e \ 



2|s(z)|e/W/ e 

= exponentially small in e, uniformly for zp — 77 < z < zp . 
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Since f(z) changes sign at Zp, we have e s {*+(^)}/ e = e f( z )/ e for z > zp, and therefore 



221 J (l + e 2 /W/ e )|s(z)| 



(184) 



= 1 



(|s(z)e/M/ £ - l)(|s(z)e^ z )/ £ e- 2 ^W/ £ - 1) 



|s(z)|e/( z )/ £ (l + e- 2 /MA) 
1 + exponentially small in e, uniformly for zp < z < zp + rj, 



where we have again used (47) and the fact that e 2 ^( z )/ c < 1 for z > zp. Similarly 

|| s ( z )|2 e 2/( Z )/e _ X | 



\B u (z)\ = \B 21 (z)\ = 



(185) 



< 



2|s(z)|e/( 2 )/ £ cosh(/(z)/e) 

|| s (^)|2 e 2/( Z )A_ 1 | 



2|s(z)|e/WA 

= exponentially small in e, uniformly for zp < z < zp + ?/. 



Therefore ( 170 1 and ( 173 1 can be combined to read 

(186) P+( z ) = « cr3 P-(z)i _<T3 (I + uniformly exponentially small in e), \z — zp\ < rj. 

Now we describe the jump conditions for P(z) across the various non-real contours pictured in Figure [8] 
(it suffices to consider those in the open upper half-plane only due to the symmetry 



(187) 



P(z) = a 1 P(z*)*a 1 



inherited from ( |70[ ) and ( 164) via the fact that g{z*) — g{z)*). First consider the arc B separating the regions 
£ and TZ in Figure [8j oriented from z = q toward z = zp. Then recalling the function h(z) defined in terms 
of g(z) by (84), we have 



(188) P + (z) = P_(z) 



e i(h+(z)+h-(z))/t 
_ e -i(h + (z)+h-(z))/e q 



P_(*) 



_ p — is/e 



z G B, 



where we have used (148). Next, let B + denote the oriented arc in Figure [8] from z = q to z — zp + r\. Then 



(189) 



P+(z)=P_(z) 



1 

e -2ifc(z)/6 y 



E(z)- a ^ 2 , zeB^ 



Similarly, let B~ denote the oriented arc in Figure [8] from z — q to z — zp — 77. Then 



(190) 



P + (z)=P_(z)£(z) 



0-3/2 



1 

e -2i7i(z)/e ^ 



z g B-. 



Finally, consider the unbounded arc of E in the upper half-plane, which we denote by C. We assume that C 
is oriented from z = —00 toward z — q. Then 



(191) 



P+(z) =P_(z)^(z) CT3 / 2 



'j e 2ih(z)/<L 

1 



E(z)- a3 / 2 , zeC. 



It follows from the inequalities (146)-(147) and the asymptotic estimate (52) that all three of the jump 
conditions (189)-(191| are of the form P+(z) = P_(z)(I + 0(e)) uniformly for z on the relevant contours 
bounded away from the common endpoint z = q. In fact, for such z the dominant contribution to the error 
comes from the matrix factors E(z) ±a --l 2 as e ±2lh ^/ e is in each case exponentially small. Moreover, due 
to the estimate (149), along the unbounded contour C, we have that P_(z)~ 1 P+(z) — I is exponentially 
decaying as z —¥ 00 for each sufficiently small e > 0. 
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5.2. Parametrix construction. Let D denote a small open disk of radius independent of e centered at the 
point z = q. We now define an ad-hoc approximation to P(z), the global parametrix P(z): 

fP in (z), zGL>\£ 
(192) P(z) := i (7iP in (z*)Vi, zeJ)*\S 

[p out (z), 26C\(EUBUD*). 

We call P out (z) the outer parametrix and we call P m (z) the inner parametrix. We emphasize that there is 
no special parametrix needed in any neighborhood of the point z — zp from which the phantom poles of s(z) 
emerge into the upper half-plane. 

5.2.1. The outer parametrix. The outer parametrix is easy to write down explicitly: 



(193) 



\z) := e lR ^l {2 ^\](i{i{-z) 1 ' 2 Y 3 \J^ e- lKa3 ^ 2t \ U 



1 

71 



g-iTT/4 

e i7r/4 



u 



Ut. 



Here, i(-z) 1 ' 2 G C+ and j3(k) is the function whose fourth power is given explicitly by 



(194) 



(3(k) 4 ■■= 



[k-q 1 / 2 )^- 



el/2 



(fc-g* 1 /2)(fc + gl/2) 

where the principal branch of the square root is meant in all cases, whose branch cuts lie along the arcs of 
±(£> U B*) 1 / 2 , and for which the specific branch of the one-fourth power is selected so that /3(0) = 1. It is 
easy to see that /3(k) is analytic along the imaginary £;-axis, and since /3(fc) 4 > for imaginary k, (3(k) > 
for imaginary k as well, leading to the conclusion that /3(oo) = 1. This allows us to asymptotically expand 
j3(k) for large k, leading to the corresponding expansion of P out (z): 



(195) 



'(*)=! + 



S{qrV2} e V< 

-9f{gV2} e -i«/« o 



i(- z y/ 2 



+ 0(z- 1 ), z 



It is easy to check that P out (z) is analytic for z G C\(BL)B* UR+), and that it is bounded independently of 
e uniformly for z € C\ (DUD*) and has determinant one. Also, P out (0) = I, p out (z*)* = a 1 P out (z)cr 1 , and 

e lK / £ 



(196) 



■pout 



(z) = P°»\z) 



z £ B, 



which should be compared with (188), while P^ ut (z) = i CT3 P° ut (z)i ff3 for z > 0, which should be compared 
with |178| , ([179]), and ( p6| ). 

5.2.2. T/ie inner parametrix. We now construct the inner parametrix in terms of Airy functions [1 , following 
closely the discussion in section 5.2 of [2]. Note that as a consequence of the continuity of h(z) at z = q and 
the condition (148) we have 2h(q) = K. Then, using (150) shows that the equation 

(197) y 3 = (2ih(z) - in) 2 

defines three different univalent functions y of z in a neighborhood of z = q each of which maps z — q to 
y = and which differ by factors of the cube roots of unity; we choose the branch y = y(z) for which the 
image of B D D is a segment of the negative real y-axis that abuts the origin. We then set C : = e ~ 2 ^ 3 y( z )- 
Then for z G B + U 0_ U C we have 2ih(z) — in = — e£ 3 / 2 where the principal branch of the 3/2 power is 
meant. Note that in terms of £ we may write the outer parametrix in the form 

(198) p°"t( z ) =e iW(2e) H ^) e ^/6^3/4 U t e -i«<73/(2£) ! 

where H(z) is an e-independent analytic matrix function with determinant one, defined for z G D by 

(199) H(z) :=U[y(z)- 1 / 4 / 3(i(-z) 1 / 2 )] CT3 , z G 

That H(z) is analytic follows from the fact that y(z)~ 1 / i j3(i(— z) 1 / 2 ) extends from 13 \ B to all of D as a 
single- valued non- vanishing analytic function. Now for convenience write £ := (|) 2 ^ 3 C arL d define a matrix 
in the £-plane by the formulae: 



(200) 



Z(C) := V2tt 



o- 3 /6 



- 3 ™/ 4 Ai'(£) e ll7 "/ 12 Ai'(£e- 2 ™/ 3 ) 



*/ 4 Ai(£) 



i/12 



Ai(£e 



-2iri/3 



e 2?3/2CT3/3 , < arg(C) < 



2?r 

T' 



3.1 



(201) Z(C) := V2tt 



0-3/6 



e -57ri/12_^//^ e 27ri/3j e ll7ri/12^^ e -27ri/3V 
e -7«/12 Ai ^ e 2«/3) e «/12 Ai ^ e -2«/3) 



27T 



2^V 3 / 3) ^ <arg(C)<7r) 



(202) Z(C) := V27r~ 



<t 3 /6 



e ll7ri/12 A j'^ e -27ri/3^ e 7iri/12^ f^m/Z\ 
e ,ri/12 Ai ^ e -27ri/3) 



2 « 3/2 ^/ 3 , -7T < arg(C) < 



2?r 

'T' 



(203) Z(C) := V2vr ( - 



4 \ -3/6 r e -3™/4 Ai ^) e 7 7 rVl2 Ai '^ e 2 3 ri/3y 



e -7rV4 Ai ^) e 57ri/12 Ai ^ e 2,ri/3) 

From the well-known asymptotic behavior [T] of Ai(£) and Ai'(£), it follows that 

o(r 3 / 2 ) o(r x ) 



o C 3/2„ /o 27T , . 

e 2 « ^/ 3 , -—< arg(C) < 0. 



(204) 



z(c)uc CT3/4 = 1 + 



o(r 2 ) o(r 3/2 )j 

as £ — > 00 in all directions of the complex plane. We define the inner parametrix by setting 
(205) P in (z) := e iK<T3/(2e) H(z)e CT3 / 6 Z(Oe-'" tCT3 /( 2e )D(z), z € D \ E, 

where D(z) is a near- identity diagonal matrix factor given by 



(206) 



D(z) 



z e Dn(CuK) 



E{z)~ a3 ^, elsewhere in D \ E. 
It is a calculation to confirm that the following jump conditions hold, which should be compared with 



( 188 1-( 191 1 



(207) 



P%{z) = PL n (z) 



e lK / f 

_ e -iK/e o 



(208) 



P?(z) = P!_ n (z) 



1 

e -2ih(z)/e -y 



(209) 

and 

(210) 



P' + n (z) = P™(z)E(z) 



P?(z) = P^z)^)" 3 / 2 



03/2 



1 

e -2ih(z)/e 1 



1 e 2ih(z)/t 

1 



^(z)-' 73 / 2 , zeCni). 



That is, the inner parametrix satisfies exactly the same jump conditions satisfied by P(z) in a neighborhood 
D of z = q. 

5.2.3. Mismatch of the global parametrix across dD. It is not difficult to see that the inner and outer para- 
metrices match each other quite well on the boundary of the disk D. Indeed, we have 

(211) P in {z)P out (z)- 1 = e lKCT3 ^ 2£ )H(z)e' T3 / 6 Z(0UC o ' 3/4 e~ o ' 3/6 H(z)- 1 e- 4 ' t ' T3/£ • P out (z)D(z)P out (z)- 1 , 

z e dD. 

But, since D(z) = I + O(e) while P out (z) and its inverse are uniformly bounded for z £ dD, the latter 
factor is obviously I + O(e). Also, since z € dD corresponds to = 0(e 2 / 3 ), from (204) we have that 

e ^/6 Z ^JJ^ a /4 e -a a /6 = J + ( £ ) for z g Q D It foH ows that 



(212) 



pfzjr'fz)- 1 



O(e), uniformly for z £ dD. 
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Figure 9. The jump contour Ee for E(z). The disk boundaries are taken to be oriented 
in the clockwise direction, and all other contour arcs are oriented left-to-right. 



5.3. Error analysis. Let E(z) := P(z)P(z) 1 denote the mismatch between the (not explicitly known) 
matrix function P(z) and its explicit global parametrix P(z). Because P(z) and P(z) satisfy exactly the 
same jump conditions on the contours B, B* , and E n (D U D*), E(z) may be extended analytically to the 
latter contours. On the other hand, E(z) has jumps across the disk boundaries 3D and 3D* where P(z) 
has no jump. That is, E(z) may be regarded as an analytic matrix function of z E C \ Ee, where Ee is 
the contour shown in Figure [9j By combining the jump conditions satisfied by P(z) with those known to be 
satisfied by the parametrix P(z), it is a simple calculation to see that on the complex arcs of Ee as well as 
for z < 0, we have that E_(z)" 1 E + (z) — I is 0(e) in both the L°°(E E \ R+) and the L 2 (E E \ K+) senses. 
For z > we have that [z cr3 E_(z)i _<T3 ]~ 1 E + (z) —I is 0(e) (in fact, exponentially small) in both the L°°(R + ) 
and L 2 {R+) senses. Because both P(0) = I and also P(0) = P out (0) = I, we have E(0) = I. 

The unusual nature of the jump condition on M. + can be understood by introducing the matrix function 
F(w) defined by 



(213) F(to) 



i CT3 E(w- 2 )r CT3 , > 

E(w~ 2 ), < 0. 



The jump contour Ep for F(w) consists of reciprocals of the positive and negative square roots of the complex 
arcs of Ee along with the real and imaginary w-axes. It follows directly that ~V(w) — I = 0(e) in both the 
£°°(Ef) and L 2 (Ep) senses, where ~V(w) :— F_(w) _1 F + (w) for w £ Ef- Also we unambiguously calculate 
the value F(oo) = I. This problem is solved by writing F(w) in terms of a Cauchy integral 

(214) FH = I + 1/ (I + X( S ))(V( 5 )-I) rfS) weC ^ F) 



2wi 7e f s — w 

where X € i 2 (Ep) solves the singular integral equation 

(215) XW-1/ Ms)(V(s)-I) ds= 1 f VW-I w g Ep- 

Here the the notation w~ indicates that the integrals are to be evaluated for w $jL Sf and the limit of w 
approaching a point on Ef is to be taken from the right side according to the orientation of each arc. By 
standard theory [23], this equation has a unique solution whose norm is O(e). Now since *V(w) — 1 vanishes 
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to all orders as w — > along Ep, F(w) has a nonconvergent but asymptotic power series representation 

oo 

(216) F(w) -I - ^F„w n , w^O, 

with coefficients given by 



n=0 



(217) 



F — 



l 

2?ri 



(I + X( S ))(V( S )-I),s- 



Now in addition to the estimates II XI 



l2(Sf) - 0(e), ||V - I|| L oc (Sf) = 0(e), and || V - I||l 2 (s f ) = 0(e), we 
have from |l77| ) and ( fl79| ) that V(w) -I = 0( e - M /( e H 2 )) holds for sufficiently small (independent of e) |u>|, 
w € Ef- Assuming without loss of generality that |w| < 1, this implies that (V(u>) — I)ui - " -1 = 0(e n+1 ) 
uniformly for small w. These estimates combined with the Cauchy-Schwarz inequality yield that F„ = O(e) 
for all n > 0. Note that Fo is necessarily a diagonal matrix (for consistency with the jump conditions of 
E(z) for large z). 

From this information, we can now obtain the asymptotic expansion of P(z) = E(z)P(z) as z — > oo: 
taking Q{w} < we have w = -i(-z)^ 1 ^ 2 , and so E(z) = F(-i(-z) -1 / 2 ). Therefore, since P(z) = P out (z) 
for large z, using (195) gives 

(218) P(z) =I + F + f(I + F ) 



1 ds, n>0. 

r) =0(e), and || V - %2 (Sl 



IK,/ ' £ 



y{g 1/2 }e iK/e 




— Fi 



1 



i(-«)Va 





y{g 1/2 }e 

Finally, since for large z g AU A*, we have M(i(~z)^ 2 ) = N(z) = O(z) = V{z)e- l 3 { - z ^l 
that the solution of the MNLS Cauchy problem is 
(219) 



0{z 



lim 

z^AuA* 



2i(-z) 1 / 2 M^i^z) 1 ' 2 ) 
a M 22 (i(-z)V2) 



lim 

z— J-oo 
z£AUA : 



2»(-z) 1 ^ 2 Pi 2 (z) 
a F22 (z) 



l/2xJ K /e 



a 



we see from ( 67 ) 



0(e), 



This asymptotic formula has the form written in the statement of Theorem [T] where 



(220) 



A(x,t) := -^{q( Xl t) 1/2 } and S(x,t) := n(x,t). 



That the derived fields p{x,t) := A(x,t) 2 and u{x,t) := S x (x,t) solve the dispersionless MNLS system ^ 
has already been shown, since according to (155) we can express u(x,t) explicitly in terms of q(x,t), and 
hence 

(221) p(x,t) = 4j (^{qix^) 1 / 2 }) 2 and u(x, t) = - (1 - 4|g(ar, t)\) . 



We have already seen (see the discussion at the end of j |4.3| ) that these substitutions reduce the dispersionless 
MNLS system to diagonal (Riemann invariant) form, reproducing exactly the partial differential equations 
d 1 33[ ) that were proved in Proposition [5] to be satisfied by q{x,t) and q(x,t)*. Next, we show that A(x : 0) = 
A (x) and that S(x, 0) = S (x). To see this, we note that according to Proposition [2] we have q(x, 0) = %(x) 
for x < x c , where z — $(x) is the solution in the upper half-plane of the turning point equation ( |55[ ): 
this implies that the quadratic x(x; z) is a constant multiple of the factored form (z ~ i{x)){z — l{x)*) = 
(z — q(x, 0))(z — q(x, 0)*), which in turn implies the identities: 



(222) 



a 2 p (x) + -(au (x) - 1) = - [q(x, 0) + q{x, 0)* 



and 



— (au (x) 



\q(x,0)f 



Under the condition that u (x) < I /a, it is easy to see that (221) and (222) imply that 
(223) p{x, 0) = po(x) and u(x, 0) = uq(x). 

But because po(x) > 0, the condition uq(x) < 1/a follows from the inequality Q < that holds for x < x c at 
t = 0. It is obvious from (220) that A(x,t) > for all x < x c {t) and t > 0, and hence we can take a positive 
square root to find that A{x, 0) = ^lo^) — sech(a;). Since we already have u(x, 0) = uo(x), to prove that 
S^cc, 0) = Sq{x) we simply need to establish the latter identity at some point x < x c . In fact, we examine 
the limit x j" x c ; since according to ( 11 ) we have So(x c ) = 5x c + p log(cosh(a; c )), we simply compare with the 

: k(x, 0) as x f x c as given by ( 160 1; thus the problem is 
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explicit expression for the limiting value of S(x,0) 



reduced to checking the equality of two explicit functions of the three real parameters a, 5, and (x over the 
region constrained by the three inequalities ( 15 )— ( IT I . We have confirmed this equality numerically, and we 
therefore conclude that indeed S(x,0) = So(x) = Sx + fi log(cosh(a;)) holds for all x < x c . 
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